NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 



REPORT No. 382 

ELASTIC INSTABILITY OF MEMBERS 
HAVING SECTIONS COMMON IN AIRCRAFT 
CONSTRUCTION 

By GEORGE W. TRAYER and H. W. MARCH 




REPRODUCED BY 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

Ui. DEPARTMENT OF COMMERCE 
SPRINGFIELD. VA. 22161 



1931 



AERONAUTICAL SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 





Symbol 


Metric 


English 


Unit 


Symbol 


Unit 


Symbol 


Lengths 


I 

t 

F 


meter 


m 
s 

kg 


foot (or mile) __ 


ft. (or mi.) 
sec. (or hr.) 
lb. 


Time 




second (or hour) «_ 


Force 


weight of one kilogram 


weight of one pound 


Power 


P 


kg/m/s 




horsepower 


hp 

m. p. h. 

f. p. s. 




fkm/h 


k. p. h. 
m. p. s. 


mi./hr. 


! 


Im/s 


ft./sec. 



Wj Weight = mg 

gy Standard acceleration of gravity = 9.80665 
m/s2 = 32.1740 ft./sec.^ 

7n, Mass = — 
9 

p, Density (mass per unit volume). 
Standard density of dry air, 0.12497 (kg-m"'* b, 

s2) at 15° C. and 750 mm = 0.002378 c, 

Gb.-ft.-^ sec.2). 
Specific weight of "standard'' air, 1.2255 S' 

kg/m3 = 0.07651 Ib./ft.^ m, 



2. GENERAL SYMBOLS, ETC. 

mk^y Moment of inertia (indicate axis of the 
radius of gyration by proper sub- 
script). 
Area. 



s, 

G, 



Wing area, etc. 

Gap. 

Span. 

Chord. 



L, 
D, 
Do 



True air speed. 

Dynamic (or impact) pressure =2 P 

Lift, absolute coefficient Csu=^ 

D 



VI 

u 



Do 



Drag, absolute coeflBcient ^^"^^ 
Profile drag, absolute coefficient ^^^0"="^ 
Du Induced drag, absolute coefficient Cm^^ 



B, 



Parasite drag, absolute coefficient ^x>p'="'^ 
Cross-wind force, absolute coefficient 



Resultant force. 

Angle of setting of wings (relative to 

thrust line). . ^""^ 

Angle of stabilizer setting (relative to 

thrust line). 7? 



Aspect ratio. 

Coefficient of viscosity. 

3. AERODYNAMICAL SYMBOLS 

Qj Resultant moment. 
0, Resultant angular velocity. 

> Reynolds Number, where / is a linear 

dimension. 

e. g., for a model airfoil 3 in. chord, 100 
mi./hr. normal pressure, at 15° C, the 
corresponding number is 234,000; 
or for a model of 10 cm chord 40 m/s, 

the corresponding number is 274,000. 
Center of pressure coefficient (ratio of 
distance of c. p, from leading edge to 
chord length). 
Angle of attack. 
Angle of downwash. 
Angle of attack, infinite aspect ratio. 
Angle of attack, induced. 
Angle of attack, absolute. 

(Measured from zero hft position.) 
Flight path angle. 
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SUMMARY 

Two fundamental problems oj elastic stability are 
discussed in this report, which was prepared by the 
Forest Products Laboratory ^ for publication by the 
National Advisory Committee j or Aeronautics. In Part 
I formulas are given for calculating the critical stress at 
which a thin, outstanding flange of a compression mem- 
ber will either wrinkle into several waves or form into a 
single half wave and twist the member about its longi- 
tudinal axis. A mathematical study of the problem, 
which together with experimental work has led to these 
formulas, is given in an appendix. Results of tests 
substantiating the recommended formulas are also pre- 
sented. In Part II the lateral buckling of beams is 
discussed. The results of a number of mathematical 
studies of this phenomenon have been published prior to 
this writing, but very little experimentally determined 
information relating to the problem has been available 
heretofore. Experimental verification of the mathemat- 
ical deductions is supplied in this report. 

INTRODUCTION 

Designing for the greatest load with a given amount 
of material in a compression member generally leads 
to the distribution of material at the greatest possible 
distance from the neutral axis of the member. The 
extent to which such distribution can be carried is 
limited by the possibility of secondary failure. Com- 
pression members with relatively wide and thin out- 
standing parts may fail through local wrinkhng or 
through twisting about the longitudinal axis at loads 
considerably less than those that would be expected 
to cause the more common failures of crushing for 
short lengths or flexure for longer lengths. When 
such a compression member does fail, a thin, outstand- 
ing element may either break up into several waves 
(wrinkle) or may buckle into a single half wave, de- 
pending upon the length and the torsional resistance 
offered by the member of which it forms a part. Such 
action has been observed for years. (References 2, 
14, 15, 18, and 21.) 

Again, the strength of a beam increases more rapidly 
with depth than with thickness, and consequently in 

» Professor of mathematics, University of Wisconsin. 

2 Maintained at Madison, Wis., in cooperation with the University of Wisconan 



aircraft, where weight is such an important matter, 
designers customarily use comparatively deep, narrow 
beams. The ratio of depth to breadth, however, has 
been kept within certain arbitrary or conventional 
limits in commercial practice, because of the well- 
known fact that a beam much deeper than it is wide 
may buckle laterally and twist before it will break 
by bending in a vertical plane. As a matter of fact, 
there is for each condition of loading and support a 
critical buckling load for such a beam just as there is 
a critical Euler load for a long column. 

Either buckling or twisting or both are Hkely to 
occur in one member or another of an aircraft struc- 
ture, and hence failure of a particular member may 
be either in a normal type of bending or compression 
resulting from the normal loading or through lateral 
buckling, wrinkhng, or twisting under stresses having 
their origin in the normal loading. Means of esti- 
mating the stress at which elastic instability is likely 
to occur have therefore become necessary in the close 
designing of the present day, in order to provide 
against secondary failure. Realizing this, the Bu- 
reau of Aeronautics, Navy Department, financed an 
investigation of fundamental phases of elastic insta- 
bility to be conducted by the Forest Products Lab- 
oratory. Wood was used in the experiments, not that 
the problem is limited to any one material, but be- 
cause of the convenience with which test specimens 
can be made of wood. 

The wrinkling and twisting problem has been investi- 
gated mathematically for homogeneous, isotropic ma- 
terials, and useful results have been obtained, notably 
' by Timoshonko. (References 17 and 21 .) This report 
reviews the general theory, adds an analysis that ap- 
plies to nonisotropic material such as wood, and dis- 
cusses the diminution of the critical stress caused by 
the elastic giving of the material at the base of the 
flange. The exact mathematical approach to the 
problem leads to rather complicated results; through 
consideration of test data, however, these results can 
be simply expressed for problems of practical interest. 

The alUed problem of the lateral stability of deep 
beams has already been investigated rather fully from 
a mathematical standpoint. The results of such work 
have been pubhshed by Michell, Prandtl, Timoshenko, 
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and others. (References 9, 11, 13, 17, 20, and 23.) 
This report adds experimental verification of the re- 
sults already obtained. 

TEST MATERIAL 

Test specimens were made of Sitka spruce cut in 
Oregon and shipped in log form to the Forest Products 
Laboratory where the wood was sawed into lumber, 
marked, and seasoned. As a result of this procedure 
the history of each piece and its location with respect 
to others in the same log were known. Part of the 
lumber was immediately kiln-dried after sawing and 
part was left to air-dry. Specimens were made from 
both the kiln-dried and the air-dried stock. 

In selecting pieces for test specimens, the usual 
Army and Navy specifications were adhered to with 
an additional limitation as to knots and pitch pockets 
in that none was permitted, no matter how small. 



The elastic properties of the material in the vario n 
planks from which the major test specimens were taken 
were deternuned by testing small control specimens 
cut from the same planks and so located as to be 
representative. In certain instances it was possible 
to accomplish the same result by cutting the control 
specimens from uninjured portions of the major test 
specimens after the main test had been completed. 
In other instances such properties as the stiffness in 
bending and the torsional rigidity of major test speci- 
mens were determined by a secondary test of the major 
specimens themselves either before or subsequent to 
the main instability test. In such secondary tests the 
stresses were kept well below the elastic limit and 
when they were made the usual control tests served 
only as a check. 
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PART I 

THE WRINKLING AND TWISTING OF COMPRESSION MEMBERS HAVING THIN, OUTSTAND- 
ING FLANGES 



METHOD OF TEST 

WRINKLING TESTS 

Two principal series of wrinkling tests were made on 
compression members having thin, outstanding flanges. 
In one series, a number of specimens, all having a 
single flange of the same size, were tested under a com- 
pressive load and the half wave length and the load at 
which w^rinkling started were recorded. The out- 
standing flanges were then reduced in width a given 
amount with the thickness left as before and the speci- 
mens were again tested. This procedure was con- 
tinued until the widths had been so reduced that 
wrinkling did not occur. 

In the other principal series of tests the width ol" 
flange was kept constant and the thickness was re- 
duced after each test. Several specimens were used in 
order to obtain reliable averages for the half w^ave 
length and the wrinkling stress corresponding to each 
thickness. Figure 1 shows a specimen in the testing 
machine. 

In addition to the two principal series of tests, a 
number of tests were made on built-up U, I, and +, 
sections under axial compression. 

TWISTING TESTS 

The set-up for the twisting tests is shown in Figure 2. 
Extension screws w^ere attached to an ordinary 4-screw 
testing machine in which specimens up to several feet 
in length could then be handled. This set-up was used 
only to obtain maximum load. To obtain a load- 
twist curve, a 2-screw machine was used, one that could 
take specimens up to about 12 feet in length without 
the use of extension screws. A pointer approximately 
3 feet in length was attached to one flange and in some 
instances to tw^o flanges. As the colunm twisted, the 
end of the pointer passed over a plane table supported 
from the base of the testhig machine and when incre- 
ments of load w^ere read by the operator at the balance 
beam the position of the pointer was marked and the 
load set opposite such marking. 

Prior to the twisting test each specimen was tested 
in torsion in order to obtain the torsional rigidity of 
the member. The stresses were kept well within the 
elastic limit during this test. 



ANALYSIS OF THE WRINKLING AND TWISTING 
PROBLEM 

The failure of compression members that contain 
wide, outstanding parts, as illustrated in Figure 3, 
may be brought about through wrinkling of the out- 
standing parts themselves instead of through the 
normal failure of the member as a w^hole, if the out- 
standing parts are suflSciently thin. When such wrink- 
ling occurs, the outstanding flange may either break 
up into a single half wave or into more than one, 
depending upon the torsional rigidity of the member 
and the fixity of the flanges. If an outstanding flange 
projects from a member that is high in torsional stift- 
ncss, wrinkling into several waves is likely to occur 
if the ratio of the outstanding width to the thickness 
of the flange is great. On the other hand, if the 
torsional stiflness is not great, the outstanding part or 
parts may form into single half waves and twist the 
member about its longitudinal axis. The critical 
values of the stresses at which one or the other type 
of buckling occurs are discussed in the following 
parasz:raplis. 

WRINKLING 

A mathematical approach to the wrinkhng problem 
is given in tlie appendix, where it is shown that the 
cntical value of the compressive stress p for a plate 
perfectly fixed along one edge, free along the opposite 
edge, and simply supported along the ends to which 
the load is applied is given by 

V = kE^^2 (1) 

in which h is the thickness of the plate, h its width, 
E the modulus of elasticity of the material, and k a 
coefficient depending upon the ratio of the length of 
plate a to the width h. 

The appendix shows further that for structural steel 
the calculated minimum value of A: is 1.16 and corre- 
sponds to a ratio of a to 6 of 1.6 or a multiple thereof. 
(Reference 21.) The theoretical formula for the mini- 
mum critical stress for steel would therefore be 

2)=1.16£p (2) 
7 
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The mathematical analysis, as already pointed out, at 
first assumes perfect fixit}^ at the base of the outstand- 
ing flange, a condition probably never realized in actual 
practice. Coupoquontly a critical stress much lower 




Figure 1 .— The wrinkling under load in the testing machine 
of a compression member having a single thin, outstanding 
flange 



than that predicted by the theory is to be expected. 
Roark, who used specimens like B and C of Figure 3, 
in which the outstanding flange was clamped between 
angles, foimd that the formula 

p = 0.6£:j; (3) 



represented his experimental results reasonably well. 
(Reforencos 14 and 15.) The reduction of the coeffi- 
cient from 1.16 to 0.6 can be attributed to the lack of 
perfect fijcity at the base of the flange. Even when an 
outstanding flange and the rigid back from which it 
projects are nil in one pioco, perfect fixity at the base 




Figure 2.— The twisting under load in the testing machine 
of a compression member having several thin, outstanding 
flanges 



of the flange can not be assumed. There is an elastic 
giving at the base of the plate and also in every device 
used in an attempt to obtain perfect fixity. Hence the 
exact coefficient that should be used for steel and other 
metals remains to be determined by experiment. A 
discussion of the situation for wood follows. 

The appendix shows that, on the basis of the differ- 
ential equation of a nonisotropic elastic plate, such as 
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wood, a critical half wave length and a critical stress 
may be calculated. The same mathematical work, 
liowever, also shows that the values of the half wave 
length and the critical stress vary over a wide range as 
the inclination of the growth rings to the faces of the 
outstanding flange varies from 0° to 90°. 

The fact that perfect fixity at the base of the flange, 
as at first assumed in the mathematical study, can not 
be obtained is true particularly of wood, which further 
complicates the problem. The stresses at the base of 
the flange resulting from the bending of the flange are 
acting perpendicularly to the grain of the wood, the 
direction in which wood is weakest. 

The appendix shows that the critical stress for a 
c^uarter-sawn flange of spruce perfectly fixed at the 
edge is 

2? = 0.228£'p 

For a similar flange with growth rings at 45° to the 
faces the critical stress is 



p = 0A17E^. 



Because of the elastic giving of the material at the 
base of the flange, however, there is a great reduction 
in the actual critical stress. Furthermore, this elastic 
giving tends to decrease the difference between the 
critical stresses for flanges with growth rings at 45° 
and 90°, respectively. Tests gave as the reduced 
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Figure 3. — Typical cross sections of compression members that have 
wide, thin, outstanding parts 

coefficient 0.07 for spruce flanges, and the expression 
for the critical wrinkling stress then becomes 



p = 0.07E 



(4) 



Probably this coefficient may be applied to other 
species without appreciable error. 

In Figure 4, in which wrinkling stress is plotted 
against the ratio of flange width to thickness, are 
shown the results of some actual tests. Each circle 
represents the average of from 4 to 18 values. The 
results have been adjusted by direct proportion to 
36276—31 2 



correspond to a modulus of elasticity along the grain 
of 1,000,000 pounds per square inch. The full line is 
the locus of the expression 

/ 2 

p = 0.07Xl,600,000p (5) 

No record of the angle between the growth rings and 
the faces of the flange, the importance of which has 
been mentioned, was made at the time of test, but 
full-section blocks from many of the test specimens 
were saved and the angle was subsequently measured. 
The direction of the rings ranged from 45° to 90°, as it 
does in what may be called commercial edge-grain 
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Ratio of flange widfh to tliickness. b/ti 

-The relation between the ratio of flange width to thickness and the 
wrinkling stress of thin, outstanding flanges 



(quarter-sawn) stock. The test specimens, therefore, 
represent what would be found in actual practice. The 
variation in the test results is accounted for by the 
variation in the direction of the growth rings and the 
difficulty of determining accurately just when wiin- 
kling started. 

Since the phenomenon of wrinkling is one to avoid in 
good design, it is unnecessary to calculate the critical 
stress with extreme precision. Merely a fair approxi- 
mation of the critical stress is sufficient to make sure 
that for the width and thickness of flange used the 
critical wainkling stress will exceed the primary stress 
expected from the normal loads. Slightly superior 
design in this regard will seldom mean an appreciable 
sacrifice in load-weight ratio. 

Length of outstanding flange. 

The coefficient k in the expression for critical wrin- 
kling stress is a minimum when the ratio of the length 
of plate a to the critical half wave length c is an inte- 
gral number. If the plate is short and ajc is not an 
integer, the critical stress may be considerably greater 
than that given by the formula because the flange can 
not then break into the ideal half wave length. If the 
length is great, that is, if ajc is greater than 2 or 3, and 
the ratio ajc is not an integer, however, the critical 
stress will be only slightly above that given by the 
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fonnula, since the plate can then break mto a half 
wave length very close to the ideal. In either case, the 
fonnula will give values on the side of safety. For 
greater detail see Tables VI, VII, VIII, XIII, and XIV 
and Figures 23 and 24 in the appendix. 



TWISTING 



It is shown in the appendix that the critical buckling 
stress for a long steel plate simply supported along one 
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12 16 20 0 

Angle of twist at center of column, degrees 

Figure 5.— The relation between angle of twist at the center of a column and axial compressive load for 

various cruciform cross sections 

side, free along the other side, and simply supported 
at the ends, to which load is applied, is given by 



V-- 



(6) 



when Poisson's ratio is taken as 0.3. 

For spruce the coefficient of equation (6) becomes 
0.044 which, as explained in the discussion of equa- 



tions (87), (88), and (89) of the appendix, is a good, 
average figure for this species. This value of the coeffi- 
cient probably may also be applied to other species 
with sufficient accuracy. The critical stress is then 
given by 

:P = 0.044£'p (7) 

If a member with a section like D of Figure 3 is sub- 
jected to compressive stress, the outstanding flanges 
will usually form into a single half 
wave at a certain critical stress and 
in so doing will twist the menib(M" 
about its longitudinal axis. When 
such action occurs, the outstanding 
elements are essentially acting as 
plates simply supported on one side 
and free along the opposite side, and 
formulas (6) or (7) are used to calcu- 
late the critical stress. 

Members with I, H, or U sections, 
such as E and F of Figure 8, likewise 
may twist under compressive loads if 
the torsional rigidity of the section 
is not great. If the torsional rigidity 
is made large by using generous fil- 
lets or, as with a U section, by mak- 
ing the back considerably heavier 
than the legs, failure through wrin- 
kling into several waves may be 
brought about and the critical stress 
in such cases must be computed by 
the formulas applying to that phe- 
nomenon. 

Actually, the rigidity of the mem- 
ber may be such that failure will 
take place at a critical stress inter- 
mediate between the minimum twist- 
ing stress and the wrinkling stress, as 
pointed out in the appendix. It is 
extremely difficult, however, to cal- 
culate accurately the coefficient for 
the intermediate conditions. Con- 
sider for the moment wood mem- 
bers with a section like D of Figure 
3. With no fillets at the junction of 
the four legs, the coefficient 0.044 
was found to apply. As fillets were 
added, the critical stress increased in 
practically the same ratio as the tor- 
sional rigidity. A U section, such 
as F of Figure 3, will twist at a 
stress corresponding to a coefficient of 0.044 if the 
back and the legs are of the same thickness. If the 
thickness of the back is increased or if fillets are added, 
the critical stress will increase in about the same 
ratio as the torsional rigidity. 

Hence the Forest Products Laboratory recommends 
that the critical twisting stress be first calculated for 
such sections as D, E, and F of Figure 3 on the sup- 
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position that no fillets are present and that all parts 
are of the same thickness. This stress should then be 
increased by multiplying it by the ratio of the torsional 
rigidity of the actual section to the torsional rigidity 
of the assumed section. This rule apphes until the 
limiting critical stress corresponding to the coeflScient 
0.07 is reached. 

In Figure 5 are shown a number of cruciform sec- 
tions, some with and some without fillets. Wood 
columns having these sections failed through 
twisting about a longitudinal axis. Accom- 
panying each section is a graph showing the 28000 
relation between axial load and the angle of 
twist for the column corresponding to it. 
The horizontal dotted lines in these graphs 
are drawn at the critical loads calculated in 
accordance with the preceding recommenda- 
tions. 

For example, the critical *fetress for column 
T-25 (fig. 5) without fillets is given by 

^-0.044X1,901,000 [3-^^^^^^ ^ 

= 2,072 pounds per square inch. 

The area is 6.76 square inches and the criti- 
cal load becomes 

^ = 2,072 X 6.76 = 14,000 pounds. 

As a further illustration, the critical twist- 
ing stress for column T-25 (fig. 5) with 
M-incli square fillets is calculated thus: 

The torsion constant K for the section 
without fillets is 

/iL =2 X 0.318 X 6.936 X (0.506)^ = 0.572. 

For the section with fillets K must be cal- 
culated in three parts — the first part is the 
value Kx for the square central portion of 
the column section, the dimensions of which 
are 1.006 inches on each edge; the second 
part is the total value K2 for the four rec- 
tangles projecting from the square center; 
and the third part is the increase K3 caused 
by the four junctions. (Reference 22, p. 26, 
and 1929 annual report, p. 696.) The junc- 
tions are treated as T junctions and the bar 
of each T is taken as half of the square cen- 
ter. The torsion constant is then the sum 
of the parts, which are calculated as follows : 
(1.006)^ 
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Actual tests of the specimens, made prior to the 
twisting tests, yielded a ratio of 1.29. 

p = 2,072 X 1.26 = 2,610 pounds per square inch 
Area with fillets = 7.01 square inches 
Critical load P = 2,610 X 7.01 = 18,300 pounds. 

In figure 6 are shown a number of U sections of 
columns that failed through twisting about a longi- 
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Figure G. 



7.11 



= 0.144 



J?2 = 2 X 0.315 X 2 X 2.965 X (0.506)^ = 0.484 
iiCa =4 X 0.15 X (0.629)^ =0.094 

K 

Then 

Torsional rigidity with fillets 0.722 
Torsional rigidity without fillets" 0.572 



= 0.722 



= 1.26. 



4 8 12 0 4 8 12 

Angle of fwisi at center of column, degrees 

-The relation between angle of twist at the center of a column and axial compressive 
load for various channel cross sections 



tudinal axis. Accompanying graphs show the relation 
between the axial load and the angle of twist. The 
horizontal dotted lines are drawn at the critical loads 
calculated by formula (7). 

The agreement between tests results and calculated 
results as shown in Figures 5 and 6 is considered quite 
satisfactory. 
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Effect of length. 

Ill airiviiig at the coefficients 0.385 for steel and 
0.044 for spruce, which arc used in the critical-stress 
formula for free twisting, the length of plate was 
assumed as several times the outstanding width. This 
assumption gives the lower limit for the critical stress. 
As the length is decreased to less than five or six times 
the width, these coefficients increase appreciably. 
Consequently, if the legs of a channel section, for 
example, are supported at intervals as by bracing and 
the distance between points of support is less than 
five or six times the width of the legs, the actual 
critical stress will be higher than that given by the 
proposed formulas. 

CONCLUSIONS FOR PART I 

Thin, outstanding flanges of compression members 
under load may buckle into several waves or may 
buckle into a single half wave, in which event they 
will tend to twist the member about its longitudinal 
axis. 

If both the length and the torsional rigidity of the 
member are great such flanges will buckle into several 
waves (wrinkle) and the critical stress for spruce 
flanges is then given by 

If the torsional rigidity of the member is not great, 
the thin, outstanding flanges will twist the member. 



Under such rigidity the flanges may be regarded as 
plates simply supported on three edges and free along 
the fourth edge. The critical stress for such a spruce 
plate is given by 

Although the coefficients m the preceding formulas 
were obtained from the test of spruce flanges, the rela- 
tions among the elastic constants for the various 
species are such that the coefficients may be expected 
to apply to all aircraft woods with safety. 

Members having sections as shown in Figure 3 will 
twist under axial compression if the junction of the 
main elements is not strengthened with fillets. If 
generous fillets are used or if part of the main elements 
of the section are made heavier than the rest, the thin, 
outstanding elements may either wrinkle or twist the 
meml)er, this depending ufton the amount of torsional 
rigidity added. Elastic instability, therefore, may 
occur at a stress intermediate between the critical 
stresses corresponding to the coefficients 0.044 and 
0.07. Intermediate critical stresses may be calculated 
by the rules given in this report. 

Failure through local buckling can occur only when 
the critical stress is less than the stress required to cause 
primary failure. 

Further conclusions, including calculated coefficients 
for steel, follow the mathematical appendix. 
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ELASTIC INSTABILITY OF MEMBERS HAVING SECTIONS COMMON IN AIRCRAFT 

CONSTRUCTION 



PART II 

THE LATERAL BUCKLING OF DEEP BEAMS 



METHOD OF TEST 
VARIATION OF FACTORS AFFECTING THE BUCKLING LOAD 

In order to determine to what degree certain factors 
affect the critical load for lateral elastic insta})ility of 



The loading device consisted of live parts. A rod 
with an upset central portion passed through the beam 
at the neutral axis. The upset portion was threaded 
at each end so that tlie rod could be centered in the 




Figure 7— The set-up of the test for lateral elastic 

deep beams, tests were made in which all factors except 
one were held constant while the isolated factor was 
varied. In these tests the beams rested on two sup- 
ports with their ends held vertical and clamped against 
lateral rotation but free to rotate in a longitudinal- 
vertical plane as the beam deflected. Load was 
applied at the center by means of the rod-and-bar 
framework shown in Figure 7. 



instability of a single beam under center loading 

beam by means of two nuts, which were drawn snug 
against the sides of the beam during test. Slotted 
bars, the lateral positions of which were fixed by V's 
in the upset rod, connected each end of the rod to the 
ends of an evener bar and from the center of this 
evener bar a tiebar passed through the movable head 
of the testing machine and was pin-connected to it 
on the under side. All connections other than the 
pin connection mentioned were knife-edge. 

13 
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^BALL BEARING 



BALL BEARING- 



FRCCLY TWISriNQ 
AND SrVINQIN^ AfTfrtS 




FREELY TVS/ 1 ST I NQ 
AND SWINGIN<} ARMS 



BASE or MACHtNC 



Figure 8.— The set-up of the test for critical buckling load under constant bending moment with 
the ends of the single beam held vertical and not restrained laterally 



When load was applied by 
lowering the movable head, the 
beam could buckle freely to one 
side or the other. The set-up was 
not considered satisfactory until 
tTie beam budded to one side and 
then to the other with the slightest 
adjustment of the rod by means 
of the two nuts. 

TESTS OF SINGLE BEAMS UNDER VARI- 
OUS LOADING CONDITIONS 

Three different loading and fix- 
ity conditions were chosen to dem- 
onstrate the applicability of the 
formulas recommended for the 
calculation of critical buckling 
loads. These conditions were: 
First, constant bending moment 
with the ends of the beam held 
vertical and not restrained later- 
ally; second, constant bending 
moment with the ends of the beam 
hold vertical and restrained later- 
ally; and third, a concentrated 
load at the center of a beam that 
rested on two supports with its 
ends both held vertical and re- 
strained laterally. 

Constant bending moment with- 
out lateral fixity was obtained 
by considering only the portion 
of a beam that was between two 
symmetrical loads. A total span 
of 14 feet was used and the two 
symmetrical load points were 60 
inches apart. In order to pei'mit 
the beam to swing freely, both 
supports and loads were applied 
through members, 16 feet long, 
that were free to swing and twist. 
The beams were wedged into these 
long members, which were slotted 
and of sufficient rigidity to hold 
the beams vertical. The two 
loading members were attached 
to an evener timber, which in 
turn was attached to the mov- 
able head of a testing machine 
with a tie bar. The set-up re- 
quired head room of approxi- 
mately 35 feet. A diagrammatic 
sketch of this set-up is shown 
in Figure 8. 

Constant bending moment with 
lateral fixity was obtained by using 
a symmetrical 2-point loading and 
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again considoring only the portion of the beam between 
the loads. For this condition, improvised extension 
wings were put on a 30,000-pound capacity testing 
machine that permitted spans up to 16 feet. Load 
was applied at tw^o symmetrical points, in some tests 
5 and in other tests 6 feet apart. In order to obtain as 
complete lateral fixity as possible at the load points, 
lateral, horizontal, pin-connected tie rods were attached 
to the beam at intervals between the load points and 
the supports. In addition, pieces 1% inches thick and 
about 6 inches deep were clamped to both sides of the 
beam from each load point outward and well toward 
the support. Figure 9 shows this assembly. 



at the supports but because of resting on ball bearings 
were not restrained laterally. Figure 10 shows a 
panel before test. 

ANALYSIS OF THE LATERAL BUCKLING PROBLEM 

A mathematical analysis of the lateral elastic insta- 
bility of deep rectangular beams leads to the following 
general expression : 

TP Itpt mz 

(8) 



in which 

P = the critical buckling load 

£'=the modulus of elasticity along the grain 




Figure 9.— The set-up of the test for critical buckUng load under constant bending 

The third method of test, namely, the application 
of a concentrated load at the center of a beam resting 
on two supports with its ends held vertical, was 
identical with the test procedure described under the 
lu-ading. Variation of Factors Affecting the Buckling 
Load. 

TESTS OF PANELS 

Panels consisting of two beams held together with 
ribs were tested in two ways. The first method was to 
suspend the two beams on hanging supports 16 feet 
long and to apply load to each beam at two symmetri- 
cal points as just described for the testing of single 
beams under constant bending moment with ends 
held vertical and not restrained laterally. The second 
method was to support the two beams on four ball 
bearings and to apply a imiforinly distributed load 
over the ribs themselves; in doing this strips were 
laid on the ribs upon which cans filled with sand were 
placed. The ends of the beams were held vertical 



moment with the ends of the single beam held vertical and restrained laterally 

/2 = the moment of inertia about the principal vertical 



axis 

G=i\\Q moduhis of rigidity in torsion 
K=i\\e torsion constant of the section 
L = the span 

F=a. constant depending upon the loading and fixity 

conditions. 
(References 9, 11, 13, 17, 20, and 23.) 

If b is taken as the width of beam and d the depth, 
I2 in equation (8) becomes 

and the torsion constant K is expressed as follows: 

K=pdb' (9) 

in which is a constant depending upon the ratio of 
d to b. Table I gives the values of 13 for various 
ratios of d to b. 
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TABLE I 

THE FACTOR (i FOR CALCULATING THE TORSIONAL 
RIGIDITY OF RECTANGULAR PRISMS 



Ratio of 


P 


Ratio of 


o 
P 


sides, d/b 


sides, d/b 


1.00 


0. 14058 


2.25 


0. 24012 


1.05 


. 14744 


2.50 


. 24936 


1. 10 


. 15398 


2. 75 


. 25696 


1. 15 


. 16021 


3. 00 


. 26332 


1.20 


. 16612 


3. 50 


. 27331 


1. 25 


.17173 


4.00 


.28081 


1.30 


. 17707 


4. 50 


. 28665 


1.35 


. 18211 


5.00 


. 29135 


1.40 


. 18690 


6. 00 


. 29832 


1.45 


. 19145 


7. 00 


. 30332 


1.50 


. 19576 


8. 00 


. 30707 


l.fiO 


. 20374 


9.00 


. 30999 


1.70 


. 21093 


10. 00 


. 31232 


1. 75 


.21428 


20. 00 


.32283 


1.80 


. 21743 


50. 00 


. 32913 


1.90 


. 22332 


100.00 


.33123 


2.00 


.22868 


00 


.33333 



Figure 12 shows the results of one representative series 
of these tests. The circles represent test vahies and 
the full line is the locus of equation (11). Again the 
agreement between actual test results and theory is 
considered good. 

In the third series of tests, the span L was varied 
while all other factors were held constant. The buck- 
ling load for this condition reduced to 



p=9 



(12) 



In Figure 13 are shown the results of two representa- 
tive series of these tests. Again the circles represent 
actual test values and the full lines the respective loci 
of equation (12) for the two ])eams selected. 




Figure 10.— A ribbed pauel before test 



In the first series of tests to check the relation of the 
various factors in the general equation, all factors except 
the depth of beam {d) were held constant. The buck- 
I ng load then reduces to 



(10) 



in which C\ is a constant. In Figure 11 are plotted the 
results of four series of tests in which d was varied while 
all other factors were held constant. The ch'cles repre- 
sent the actual loads and the full lines are loci of equa- 
tion (10). The agreement is considered very satis- 
factory. 

In the second series of tests, the widtli 6 was varied 
while all other factors were kept constant. The buck- 



ling load in this case becomes 



(11) 



The effect of the nioduhis of elasticity in bending 
could not be separated fi-om that of the modulus of 
rigidity in torsion for the purpose of checking further 
the fundamental expression, because when one is 
changed the other changes with it, and therefore 
neither could be isolated. Moreover, it was impossible 
to ascertain experimentally with wood alone the impor- 
tance of their combined effect on buckling load because 
the range over which their product varies is too limited. 
For steel, the modulus of rigidity in torsion is commonly 
taken as two-fifths of the modulus of elasticity in 
bending while for spruce it is in the neighborhood of 
one-fifteenth or one-sixteenth. Since some previous 
tests of steel beams have shown excellent agreement 
with critical values calculated by the formulas, it 
therefore appeared logical to assume that, if tests of 
wooden beams also checked values given by the formu- 
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las, the moduli of elasticity in bending and of rigidity 
in torsion are in their right relation in the formula. 
(References 6 and 9.) 

Following are formulas that apply to rectangular 
beams under various loading and fixity conditions. 
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Case 1. — A thin, deep, rectangular beam under con- 
stant bending moment M, with its ends 
not restrained. 
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Figure 11.— The relation between the lateral buckling 
load and the depth^ of beam modified by a torsion cor- 
rection factor (dV/^)» for deep, rectangular beams 

In all cases the ends of the beam are assumed to be 
vertical. An end not restrained, in the terminology 




20 40 60 

Span, inches 

Figure 13.— The relation between the lateral buckling load and 
the span, for deep, rectangular beams 

Case 2. — The same as case 1 except that the ends are 
restrained. 



_ 2t'^EI2GK 
L 



A . No lateral restraint at ends 



Plan of undeflected beam 




Plan of longitudinal axis after deflection 



Q Ends restrained laterally 




Plan of longitudinal axis after deflection 



C Lateral restraint^ at center of span 



08 



24 



.16 

Figure 12.— The relation between the lateral buckling load and 
the cube of the width of beam modified by a torsion correction 
factor ib'^-y/S), for deep, rectangular beams 

used, is held vertical but is not otherwise constrained, 
and an end restrained is both held vertical and clamped 
against lateral rotation. Figure 14 shows the lateral 
deflection of the longitudinal axis for three principal 
conditions of restraint. 

36276—31 3 




Plan of longitudinal axis after deflection 



Figure 14.— The lateral deflection of the longitudinal axis of a single rectangu- 
lar beam when the bending in a vertical plane becomes unstable and sidewise 
buckhng occurs 

Case 3. — A thin, deep, rectangular cantilever with a 
concentrated load P at the end. 



P = 



Case 4. — A thin, deep, rectangular cantilever with a 
uniformly distributed load W. . 



TF= f-, 
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Case 5. — A thin, deep, rectangular beam supported 
at the ends and carrying a concentrated 
load P dA the middle, with its ends not 
restrained. 



P = 



Case 6. — The same as case 5 except that the ends are 
restrained. 



p 2b,^'ylEhGK 
^~ U 

Case 7. — A thin, deep, rectangular beam supported at 
its ends and carrying a uniformly dis- 
tributed load W with its ends not re- 
strained. 



Case 8. — The came as case 7 except that the ends are 
restrained. 



W- ^2 

Case 9. — A thin, deep, rectangular beam subjected to 
a constant bending moment M and an 
axial thrust P\ with its ends not restrained. 



T^El2GK 



EI2 



Case 10. — The same as case 9 except that the ends are 
restrained. 



^_2jvjlEhGK L P'U 

Case 11. — A thin, deep, rectangular beam supported 
at its ends and carrying both a uniformly 
distributed load W and a concentrated 
load P at the middle, with its ends not 
restrained. 

16.9 "^28:3 "^^^2^^ 

Combinations of the preceding cases may be similarly 
expressed . 

Case 12. — A thin, deep, rectangular beam supported 
at its ends and carrying a concentrated 
load P at its middle,, with lateral support 
as by tie-rods, at the middle, and the 
ends not restrained. Such a beam 
buckles laterally in two half waves. 
(Fig. 14, C.) 



P = 



U.5^|EI2GK 



BUCKUNG FORMULAS FOR I BEAMS 

The preceding formulas require modification when 
the beam has flanges, since the lateral flexure of the 
flanges then becomes important. Following are some 
of the results obtained by Timoshenko. (References 
7, 16, 17, 18, and 20.) Two more symbols are intro- 
duced. Let 

/3 = the moment of inertia of one flange about the 

principal vertical axis 
and let 

2^ EI,h' 



2 GKL ' 



Case 13. — An I beam subjected to a constant bending 
moment M, with its ends not restrained. 



Case 14. — The same as case 13 except that the ends 
are restrained. 



Case 15. — The same as case 13 with the addition of an 
axial thrust P\ 



1/ 

EL 



Case 16. — The same as case 15 except that the ends 
are restrained. 



M=-^L Vl + 4.V./1-^ 



EI2 

Case 17. — A cantilever I beam with a concentrated 
load P at the free end. 



F^|EI2GK 



in which values of F for reciprocal values of ci\ are: 

^/ 0. 1 1 2 4 8 12 16 24 32 40 « 
F: 44. 3 15. 7 12. 2 9. 8 8. 0 7. 2 6. 7 6. 2 5. 9 5. 6 4. 0 

Case 18. — An I beam supported at its ends and carry- 
ing a uniformly distributed load ir, with 
the ends not restrained. 



^.,_F^lEhGK 

in which values of Fy for reciprocal values of and for 
three differant placements of the load, are: 



— : 0.4 4 



16 32 48 64 



160 320 



(1) F; 143.2 53.0 42.6 36.3 32.6 31.5 30.5 29.8 29.2 28. 6 28. 3 

(2) F: 92.8 36.3 30.4 27. 4 26. 2 26. 2 25. 8 26. 0 26.2 26. 5 28. 3 

(3) 221.6 78.2 59.4 43.1 40.7 38.1 36.0 34.4 32.6 31.0 28.3 

The placements of the load on the beam, numbered to 
correspond with the values of F, are: 

(1) Along the neutral axis. 

(2) On the top. 

(3) At the bottom. 



\ 
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Case 19. — The same as case 18 except that the ends 
are restrained. 

F^ EhGK 

p — 

0.4 4 8 16 32 96 128 200 400 co 

a* 

F; 488 160.8 119.2 91.2 73.0 58.0 55.8 53.4 51.2 43.3 

Case 20. — An I beam supported at its ends and carry- 
ing a concentrated load P at the middle, 
with the ends not restrained. 

^ F^EhGK 

— 2= 0.4 4 8 16 32 64 96 160 320 oo 

(1) F: 86.4 31.9 25.6 21.8 19.4 18.3 17.9 17.5 17.2 16.9 

(2) F: 51.4 20.2 17.0 15.4 14.9 14.9 15.0 15.4 15.7 16.9 

(3) F: 145. 6 50. 0 . 38. 2 30. 5 25. 5 22. 4 21. 2 20. 0 18. 9 16. 9 

As in case 18 the load is applied: 

(1) Along the neutral axis. 

(2) On the top. 

(3) At the bottom. 

Case 21. — The same as case 20 except that the ends 
are restrained. 

p F-^EhGK 
^" U 

— ^: 0.4 4 8 16 32 64 96 160 320 400 oo 

"f: 268 88.8 65.5 50.2 40.2 34.2 31.8 30.0 28.5 28.2 25.9 

Case 22. — An I beam supported at its ends and carry- 
ing a concentrated load P at the middle, 
with the ends not restrained, and the 
beam laterally supported at the middle, 
as when two parallel girders have a lateral 
connection between them at the middle 
of their span. 

^~ U 

—' 0.4 4 8 16 32 96 128 200 400 oo 
a' 

F; 466 154 114 86.4 09.2 54.6 52.4 49.8 47.4 44.5 

Case 23. — An I beam supported at its ends and carry- 
ing a distributed load IT, with the ends 
not restrained, and the beam laterally 
supported at the middle of the span. 

W= JO 

0.4 4 8 16 32 96 128 200 400 « 

a* 

(1) F: 673 221 164 125.5 100.8 79.4 76.4 72.8 69.6 65.9 

(2) F: 586 194 145 112 91.2 73.7 71.5 68.9 66.8 66.9 

(3) F: 774 252 186 141 111.2 85.6 81.6 77.0 72.5 85.9 

Again the load is applied: 

(1) Along the neutral axis. 

(2) On the top. 

(3) At the bottom. 

EXPERIMENTAL VERIFICATION OF THE BUCKLING FORMULAS 

Time and funds were not available for the experi- 
mental verification of the formulas for all the loading i 



and fixity conditions listed. Over 40 I and rectangular 
beams, however, w^ere tested imder the following condi- 
tions, which represent a considerable range for the 
fixity and the loading constant F. 

Case 1. — A rectangular beam subjected to a constant 
bending moment, with its ends not 
restrained. 

Case 2. — A rectangular beam subjected to a Con- 
stant bending moment, with it sends re- 
strained. 

Case 13. — An I beam subjected to a constant bending 
moment, with its ends not restrained. 

Case 14. — An I beam subjected to a constant bending 
moment, with its ends restrained. 

Case 21. — An I beam resting on two supports, with a 
concentrated load applied at the middle 
of the span, and the ends restrained. 

The results are shown in Tables II, III, IV, and V. 
Since the exact fixity conditions assumed in the mathe- 
matical analyses are difficult of attainment, the agree- 
ment of test results with values given by the formula 
is remarkable. We consider this agreement, together 
with the agreement for a limited number of metal 
beams, concUisive proof that the formulas are appli- 
cable to beams under actual service conditions. 

A REPORTED DISAGREEMENT WITH EXPERIMENTAL RESULTS 

The only experimental record of tests with wood 
that has come to the attention of the present authors 
is an undergraduate thesis that has been published 
as National Advisory Committee for Aeronautics 
Technical Note 232, '^The Lateral Failure of Spars." 
In this note a wide difference between actual and 
theoretical results is reported, the statement being 
made that actual loads ranged from one-half to one- 
fifth the loads calculated by the formula applying 
to the test conditions. Examination of this note, 
however, leads to the conclusion that the theoretical 
formulas were not correctly applied in two respects, 
as follows: 

1. The coefficient 16.9, which the authors of the 
note used, applies only to the conditions of case 5 of 
the present report. Their loading conditions, how- 
ever, were those of case 12, which requires a coefficient 
of 44.5. In addition, the ends of the test beam were 
under light lateral restraint, which would increase 
the coefficient to about 50. 

2. It appears that they used the moment of inertia 
about the principal horizontal axis instead of that 
about the principal vertical axis. 

Onl}^ part of the test results reported could be 
checked, since in several instances the beams were 
stressed beyond the elastic limit and stress-strain 
curves with which to modify the modulus of elasticity 
were not available, yet proper work-up of their ex- 
perimental data gives results that check with pre- 
cision the theoretical results. 
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TABLE II 

THE CONSTANT BENDING MOMENT REQUIRED TO CAUSE LATERAL BUCKLING AND TWISTING OF THIN, 
DEEP, RECTANGULAR BEAMS HAVING THEIR ENDS UNRESTRAINED LATERALLY BUT HELD VERTI- 
CALLY ALTHOUGH FREE TO ROTATE IN A LONGITUDINAL-VERTICAL PLANE 



STRESSES WITHIN THE ELASTIC LIMIT 



1 


2 


3 


4 1 5 

1 


6 


7 


8 


9 


10 


11 


Beam 


Nominal dimensions (inches) 


Eli by- 


GK by- 


L 


Buckling moment by- 


Calculation from columns— 


Test 


Calculation 


Test 


Calculation 


Test 


3 and 5 


3 and 6 


4 and 6 


R-102 


H by 6- - 


398, 000 
680,000 
984, 000 
1, 744, 000 
3, 330, 000 




75, 050 
218,800 
175,900 
328,800 
547,000 


78,000 
256, 300 
161,300 
373,200 
496,000 


60 
60 
60 
60 
60 


9, 070 
20, 200 
21, 800 
39, 620 
70,700 


9,225 
21, 870 
20, 850 
42,230 
67,300 




10, 480 
18, 850 
22, 620 
43,200 
67,500 


R-106 — 


1 by 6 






R-107 

R-109 

R-111 - 


Iby6 

m by 6 

mhye - 


938,000 
1, 667, 000 
3, 638, 000 


20,380 
41,380 
70, 300 



STRESSES BEYOND THE ELASTIC LIMIT 



1 


2 


3 


4 


5 


6 7 


Beam 


Nominal dimensions (inches) 


Corrected 
E'h 


OK by test 


L 


Buckling moment by- 


Calculation 
from columns 
Sand 4 


Test 


R-110 - 


VA by 4^^... 


1,630,000 
4, 670, 000 
2, 105,000 
6, 980, 000 


389,600 
1, 059, 000 

355,000 
1,088,000 


60 
60 
60 
60 


41,730 
116, 250 

45, 250 
133,500 


42, 350 
102,800 
42, 350 
96,250 


R-112 


2 by 6 


R-113 - 


lHby4H 


R-114 - 


2 by 6. 







All calculations were made with a slide rule. 

£= modulus of elasticity as determined from control tests increased 11 per cent to correct for shear distortion. 

JS'=secaut modulus of elasticity as obtained from a stress-strain curve. 

/j= moment of inertia of a beam about its principal vertical axis. 

0= modulus of rij^idity. 

Ar= torsion constant for the section. 

L= length subjected to constant moment. / 

TABLE TIT 



THE CONSTANT BENDING MOMENT REQUIRED TO CAUSE LATERAL BUCKLING AND TWISTING OF THIN 
DEEP I BEAMS HAVING THEIR ENDS UNRESTRAINED BUT HELD VERTICALLY ALTHOUGH FREE TO 
ROTATE IN A LONGITUDINAL-VERTICAL PLANE 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Beam 


Nominal dimensions (inches) 




Eh by 
calculation 


Eh by- 


OK by 
test 


L 


Buckling moment by — 


Calculation 


Test 


Calculation 
from columns 
6 and 7 


Test 


I-IO 

I-ll 

1-12 

I-13-- — - 

1-14 - 

1-15 

1-16 

1-17 

1-18 

1-19 

1-22 


2 by 6 by H flange by ^ web 

IVi by 6 by ^2 flange by H web 

2^4 by 5M by T« flange by ^6 web 

2 by 7 by flange by M« web 

Uii by 5 by H flange by 56 web 

VAhy 6 by H flange by H web 

1 by 6 by Vi flange by web 

1 by 6 by 56 flange by H web 

2 by 7 by H flange by '^i web 

1 by 6 by Me flange by Via web 

2% by 5yi by % flange by $6 web _ 


35.40 
35.76 
26.73 
48.58 
24.60 
36.00 
35.71 
35.58 
48.58 
35.95 
27.67 


921,000 
294, 000 

2, 709, 000 
412,000 
225,000 
242,500 

80,600 
59.050 
503,000 
71,200 

3, 209, 000 


2, 140, 000 
1.006,000 
5, 710, 000 
862,000 
492,000 
536,000 
262,000 
158,000 
1,070,000 
206,000 
5,880,000 


2. 290, 000 
1,006,000 
5, 780, 000 
878,000 
529,000 
545,000 
332,000 
221,000 
1,022,000 
219,500 
5, 955, 000 


124, 100 
126, 300 
151, 500 
12, 910 
15, 380 
17, 150 
32, 250 
13, 450 
18, 400 
22,320 
150,000 


60 
60 
60 
60 
60 
60 
60 
60 
60 
60 
60 


32, 550 
19,710 
63, 100 
9, 870 

5, 770 

6, 590 
5, 750 
3, 159 

12, 070 
3, 945 
66,700 


34, 510 
21, 000 
75, 000 
13, 450 
5, 350 
8,860 
6,210 
2, 970 
15, 120 
4, 590 
65, 500 



All calculations were made with a slide rule. 
^ = height of beam. 

A'=raodulus of elasticity as determined from control tests increased 11 per cent to correct for shear distortion. 

/s= moment of inertia of 1 flange about the principal vertical axi.s of the beam. 

/2= moment of inertia of a beam about its principal vertical axis. 

G = modulus of rigidity. 

/! = torsion constant for the section. 

L= length subjected to constant moment. 
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TABLE IV 

THE CONSTANT BENDING MOMENT REQUIRED TO CAUSE LATERAL BUCKLING OF THIN, DEEP BEAMS 
HAVING THEIR ENDS RESTRAINED LATERALLY AND HELD VERTICALLY ALTHOUGH FREE TO ROTATE 
IN A LONGITUDINAL-VERTICAL PLANE 



RECTANGULAR BEAMS 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Beam 


Nominal dimensions (inches) 


Eh by- 


GK by- 


L 


Buckling moment by — 


Calculation from columns — 


Test 


Calculation 


Test 


Calculation 


Test 


3 and 5 


3 and 6 


4 and 6 


R-101 


% by 6 - 


430,200 
398,000 
984,000 




76,600 
75, 050 
175,900 




72 
60 
60 


15, 850 
18, 120 
43, 600 






14,85) 
16, 741 
36, 061 


R-102 - 


H by 6- 




78,000 
161, 300 


18, 470 
41,720 




R-107 - 


1 by 6 - 


938,000 


40,780 



I BEAMS 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Beam 


Nominal dimensions (inches). 




Eh by 
calculation 


Eh by- 


OK by 
test 


L 


Buckling moment by- 


Calculation 


Test 


Calculation 
from columns 
6 and 7 


Test 


I-10-- - 

I-11-- 

1-13 - 

1-15 

1-16 - 

1-17 - 

1-18--- 

1-19- - - - 

1-20- - 

1-21- - ---- 


2 by 6 by % flange by H web 

2 by 7 by Vie flange by Me web 

IJ^ by 6 by H flange by H web 

1 by 6 by flange by web 

1 by 6 by 96 flange by H web 

2 by 7 by H flange by ?6 web 

1 by 6 by Mr. flange by He web 

^^ by 6 by '^'{g flange by Mb web. 

by G by % (lange by H web... 


35. 40 
35.76 
48. 58 
36.00 
35.71 
35.58 
48.58 
35. 95 
35.76 
24.50 


921,000 
294,000 
412.000 
242. 500 
80,600 
59, 0:^ 
503, 000 
71,200 
241, 300 
197,000 


2, 140, 000 
1,006,000 
862,000 
536, 000 
262, 000 
158.000 
1,070. 000 
206,000 
551,000 
425, 000 


2, 290, 000 
1, 006, ()(XJ 
878, 000 
545. 000 
332, 0(X) 
221, 0(M) 
1, 022, 000 
219. 500 
547, 000 
562, 500 


124, 100 
126, 300 
12,910 
17, 150 
32. 250 
13. 450 
18. 400 
22, 320 
27, 270 
13, 650 


60 
60 
60 
60 
60 
60 
60 
60 
60 
60 


81,6r)0 
45, 100 
34. 420 
19, 700 
13, 220 

7, 780 
41,340 

9, 360 
21. 180 
15, 7'»0 


55, 900 
35, 150 
25. 380 
14, 180 
11.3.50 

7.700 
29, 430 

8, 370 
16,610 
11,480 



All calculations were made with a slide rule. 

^ = height of beam, , , ^ ^■ 

E=modulus of elasticity as determined from control tests increased 11 per cent to correct for shear distortion. 

/3 = moment of inertia of 1 flange about the principal vertical axis of the beam. 

/2 = moment of inert ia of a beam about its principal vertical axis. 

(?= modulus of rigidity. 

i\r= torsion constant for the section. 

//= length subjected to constant moment. m a -dt 

TABLE V 



THE CONCENTRATED CENTER LOAD REQUIRED TO CAUSE LATERAL BUCKLING OF THIN DEEP I BEAMS 
SUPPORTED AT EACH END WITH THE ENDS RESTRAINED LATERALLY AND HELt) VERTICALLY 
ALTHOUGH FREE TO ROTATE IN A LONGITUDINAL-VERTICAL PLANE 



Beam 



I-l---. 
B-1.. 
B-2— 
C-l.. 
C-2--. 
A-2-.. 
6-2-65 
1-4—. 
1-5.-. 
1-6---. 



Nominal dimensions (inches) 



by 6 by H flange by H web 
1 by 6 by Ir^ flange by H web... 

1 }^ by 6 by H flange by H web. 

2 by 7 by Mb flange by Mb web 

by 5 by H flange by H web. 

by 5 by H flange by H web, 
2 by 6 by H flange by lie web.. 
IK' by 5 by H flange by H web 

1 by 6 by Me flange by lie web 

2 by 6 by Me flange by lie web 



3 


4 


5 


6 


7 


8 


9 


10 


11 




E/3by 
calculation 


Eh by 
calculation 


GK by 
calculation 


L 


a' 


F 


Buckling load by- 


Calculation 


Test 


35. 05 
34. 69 
34.22 
47.61 
24. 40 
24. 21 

34. 57 
24. 21 
35.05 

35. 05 


189, 800 
74,800 
184,500 
269, 100 
137, 100 
189, 400 
339, 500 
190. 500 
71,500 
552. 000 


423,000 
236, 200 
412, 000 
562, 000 
299, 500 
414,200 
718,000 
414,800 
212,000 
1, 171,000 


15, 040 

23. 390 
13, 010 


82 
82 
82 
82 
82 
82 
82 
82 
82 
82 


0. 0328 
.0083 
.0360 


40.8 

30.6 
41.8 


482 
338 
455 


401 
405 
450 
550 
307 
400 
720 
410 
320 
O^'O 


10,080 
12, 770 
22, 960 
12, 680 
20,220 
26, 150 


.0247 
. 0268 
.0380 
.0271 
.0092 
.0550 


37.3 
38.5 
42.4 
36.9 
31. 1 
47.8 


305 
416 
810 
418 
302 
1,244 



All calculations were made with a slide rule. 

/»= height of beam. , w u j- » 

J5;= modulus of elasticity as determined from control tests increased 11 par cent to correct for shear distortion. 

/s= moment of inertia of 1 flange about the principal vertical axis of the beam. 

/2= moment of inertia of a beam about its principal vertical axis. 

0= modulus of rigidity. 

7^"= torsion constant for the section. 

L=span. 

J Eh hi 
" '"2GKL^' 

F= multiplying factor in the lateral buckling formula, dependent upon a' 



22 



REPORT NATIONAL ADVISORY COMMITl^EE FOR AERONAUTICS 



STRESSES BEYOND THE ELASTIC LIMIT 

The calculation of a critical load that produces a 
fiber stress beyond the elastic limit is possible by 
means of the preceding formulas if the modulus for 
inelastic deformation is known. Although this modu- 
lus is a variable beyond the elastic limit, it may be 
obtained from a stress-strain diagram. Figure 15 

'0000 1 I I I I I I I I I I I I I I I I I I 



in which k is a constant that need not be evaluated 
when Figure 15 is available. The modulus below the 
elastic limit will be called E in this report and that 
above will be called E\ Although both depend upon 
the slope of the line connecting the origin with the 
stress-strain curve at the particular stress in question , 
is usually spoken of as the secant modulus. 
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Figure 15.— Stress-strain curves for spruce beams. Values taken from these curves are for use in the equation: 

Modulus of elasticity =130.7 ^ train 



L6 



shows such a diagram for a spruce beam in bending. 
From it the required modulus, for a stress either below 
or above the elastic limit, may be determined by means 
of the formula : 

Modulus of elasticity = 130.7 ^ ^^train ^^^^ 



The formula proposed by Karman and advocated 
by Timoshenko for calculating 

4:EXEi 



E' = 



(14) 



{^|E^ ^E,y 

in which Ei is the tangent modulus on the compression 
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side of the beam and E is the initial modulus, can not 
be used for wood. (Reference 16.) It can not apply 
to wood because when the maximum load in bending 
is reached the stress-strain curve for the compression 
fibers has turned downward, which means that Ei has 
become negative. In fact, before the maximum load is 
reached the tangent to the stress-strain curve for the 
compression fibers has become horizontal, which means 
that the formula would give the beam no stiffness, 
whereas it actually is still resisting an increasing load. 

Whatever the method used, more than one trial will 
have to be made in the calculation of the critical stress 
because JE' is not known until the stress is known. In 
calculating critical loads by simply substituting E^ 
in the formulas that were developed on the assumption 
that the elastic limit was not passed, two further 
assumptions are made, as follows: 

1. Passing the elastic limit does not affect the 
torsion modulus G. 

2. The decrease in £'is constant along the span. 

In investigating critical loads, four rectangular ' 
beams were subjected to a constant bending moment 
that produced lateral buckling at a fiber stress beyond 
the elastic Umit. The results appear in Table II. \ 
The corrected values of £"72 given in the table were 
obtained by multiplying the secant modulus by 
the moment of inertia h of the cross section about its i 
principal vertical axis. The calculated critical bend- , 
ing moment for the first beam listed in the second part | 
of the table (R-110) is about IK per cent lower than 
the test value, while the calculated values for the 
second (R-112) and the third (R-113) beaims are 
respectively 13 and 7 per cent higher than tile test 
values. The second (R-112) and the fourth (R-114) 
beams, which were of the same size, were made from 
adjacent planks cut from the same log. Control tests 
showed the material in R-114 to be slightly superior. 
Consequently its low test bending moment is difficult 
to account for unless the beam had become slightly 
warped before test, in which event the actual stress 
at failure would be higher than the calculated stress 
and the value of E' lower than that used. 

LOAD NOT APPUED ALONG THE NEUTRAL AXIS 

The development of the buckling formulas is 
greatly simplified by the assumption that the load is 
applied along the neutral axis of the beam, and in 
aircraft work usually no material error will normally 
be introduced by assuming such an application of the 
load. In a few of the cases for w^hich formulas are 
given, coefficients are also given for load applied along 
the neutral axis, on the compression flange, and on 
the tension flange of the beam. For the development 
of the formulas for a load placed above or below the 
neutral axis, attention is again directed to the work 
of Timoshenko and to advanced texts on strength of 
materials or applied elasticity. (References 7, 12, 
and 18.) 



BUCKUNG OF BEAMS TIED TOGETHER WITH RIBS 

When two thin, deep beams are tied together with 
ribs, in addition to carrying whatever direct load is 
normally placed upon them the ribs will act to prevent 
lateral buckling of the beams. Very often, though, 
when the direct load is transferred to the beams from 
the ribs, the ribs may be laboring to sustain the load 
already upon thorn and consequently may have no 




Figure 16— The test of a panel to show that the tendency of an axially loaded 
single spar to buckle is transmitted by the ribs to an unloaded single spar 



reserve strength left for any extra load that a tendency 
of the beam to buckle would produce. 

The first panel test was made to demonstrate the 
fact that the tendency of an axially loaded spar to 
buckle is transmitted by the ribs to the unloaded spar. 
For this test there was made a panel consisting of two 
1)2 by 6 inch spars spaced 55 inches center to center, 
four compression ribs spaced 55 inches, and drag 
wires in the three bays. No ribs were put in between 
the compression ribs. Axial load was applied to but 
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one spar, which deflected alternately in and out between 
compression ribs as the beam of Figure 14, C, deflected. 
The test was stopped at a load of 12,750 pounds with 
the panel still uninjured. The deflections were in- 
creasing rapidly at that time, and apparently the load 
was very near its maximum. AuxiUary ribs were 
then put in between adjacent compression ribs, four 
in each bay. Figure 16 shows the completed panel 
ready for test. Axial load was again applied to but 
one spar. The test was stopped at a load of 29,000 
pounds, which was very near the maximum. 

Under the conditions of the second test, in which 
all ribs were in place, the two spars act as one, the 
lateral rigidity of the panel being the combined rigidi- 
ties of the two spars. Similar tests were made by the 
Engineering Division of the War Department, Air 
Service, at McCook Field with identical results. 




Figure 17.— Wing ribs for which the degree of attachment of the ribs to the beams is adjustable 

In the next panel tests the beams were subjected to 
bending, and load was applied directly to them and 
not to the ribs. Two 1 by 6 inch rectangular beams 
subjected to a constant bending moment over 60 
inches of their length were tied together with four ribs 
spaced 12 inches center to center in the 60 inch bay. 
Constant moment was applied by using the apparatus 
shown in Figure 8, except that double the number of 
support and load rods were used. The ribs that tied 
the two beams together were as shown in Figure 17; 
they were held in place simply by the friction under 
the heads of the bolts, the holes for which were slotted. 
Obviously, if the bolts were not drawn tight the beams 
could buckle very easily, while if they were drawn 
tight twisting was practically prevented. The evener 
bar was not pin-connected to the movable head in this 
test but was rigidly attached to it, so that if one beam 



stopped taking load more was thrown upon the other. 
When the panel was assembled the bolts holding the 
cleats along one beam were drawn up tightly, while 
those along the other beam were not. The beam 
supported by the less rigid cleats quit taking load at a 
moment of 35,530 inch-pounds, while the one with the 
more rigid cleats did not buckle until it was subjected 
to a moment of 53,620 inch-pounds. The results 
show what may happen when the ribs start to fail. 
Incidentally, had the beams been held so as to restrict 
bending to a vertical plane, each should have carried 
65,550 inch-pounds and had they been free to buckle 
laterally each was calculated to sustain 21,270 inch- 
pounds. 

The next panel tested was similar except that the 
ribs were glued to the flanges. Load was applied to 
the beams as before, and failure occurred when each 
beam was subjected to a moment of 
55,600 inch-pounds. The calculated 
bending moment for each with bending 
conflned to a vertical plane was 62,800 
inch-pounds. 

The third and final step was the test 
of single bays with load applied to the 
ribs alone. (Fig. 10.) The panels were 
8 feet between supports and the beams 
36 inches center to center. Seven ribs 
of the lightened plywood type, rectan- 
gular in form, extending 12% inches be- 
yond each beam and spaced 12 inches 
.'jpart, tied the two beams together. 
The ends of the beams rested on thrust 
bearings and were held vertical during 
test. Roller bearings under the ball 
bearings at one support permitted 
movement as the beam deflected. 
Thin strips 7 feet 5 inches long, notched 
at the ribs, were laid on the ribs, and 
cans filled with sand were placed on 
them. 

For this fixity and loading the beams, which were 
rectangular and by 4 inches in cross section, should 
have buckled laterally at approximately 91 pounds 
each if unsupported by the ribs. If bending had been 
confined to a vertical plane, 970 pounds should have 
been required to break each beam. The ribs when 
supported laterally should have been good for 250 to 
300 pounds. The preceding values are calculated 
ones. 

The two beams were supported at the center by a 
cross timber resting on two jack screws, with the ribs 
supported only by the beams. A load of 735 pounds 
was put on the panel and the screws lowered. The 
beams remained in a vertical plane throughout their 
length. The timber was again brought up against the 
two beams to relieve the load and more load was 
added. No buckling occurred at 1,155 pounds when 
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the screws were lowered. Again the two beams were 
supported at the center by the cross timber and more 
load was applied. A total of 1,370 pounds was sus- 
tained by the ribs with the beams still supported at 
the center. This load, however, was approaching the 
maximum for the ribs. When the screws were again 
lowered the ribs did not have sufficient additional 
strength to resist the tendency of the beams to buckle 
and they gave way. 

The two beams, which were uninjured in this test, 
were again used in a second panel. This second panel 
was like the first in every respect, but the loading was 
somewhat different. In place of the notched 7-foot 
5-inch loading strips, short smooth strips that ex- 
tended over two and three ribs alternately were used. 
Instead of having the long strips wnth their notches 
hold the tops of the ribs in line, strips }^-inch thick 
and 2 inches wide were laid flat along each side of each 
rib and tacked at the ends and center to the short 
loading strips. In this test, as in the first, the lower 
chords of the ribs were unsupported. Because the 
short loading strips permitted freer lateral play in the 
beams, this panel failed at a lower load than the first. 
A maximum load of 900 pounds was obtained, at which 
load the lower part of the ribs buckled until the ribs 
lay almost flat against the loading strips. 

In the third and final test of this series the bottoms 
as well as the tops of the ribs were held in line and the 
same beams were used again. Ten rows of l^-inch 
commercial cotton tape were run parallel to the spars 
and sewed to the ribs. Two diagonal pieces on both 
top and bottom were then sewed to the parallel rows. 
Although this taping was hardly comparable with wing 
covering, it held the ribs in line quite well. The short 
loading strips of the previous test were again used in 
addition to the tape. 

As previously stated, the lateral buckling load of 
each spar when it was unsupported w^as calculated as 
91 pounds, which is 182 pounds for the panel. The load 
required to break each one if bending had been confined 
to a vertical plane was 970 pounds or 1,940 pounds for 
the panel. Failure occurred at a total load of 1,470 
pounds, at which one beam buckled badly and col- 
lapsed. The ribs had started to buckle somewhat, 
which permitted the one beam to buckle out of a 
vertical plane. Greater strength of the ribs or greater 
torsional rigidity of the spar would have prevented this 
buckling and twisting. A box beam of the same 
strength in bending, for example, would not have 
budded at this same load. 

The nose of an airplane wing helps to hold the front 
or deeper spar in line and the \vdng covering keeps the 
ribs in line. With this support, fairly large ratios of 
depth to breadth may be used if the ribs are made with 
just a little surplus strength. 

Some years ago, after the test of a great many 
beams in connection with a study of form factors, the 



suggestion was made that the ratio of the moment of 
inertia about the principal horizontal axds to the 
moment of inertia about the principal vertical axis be 
kept low, below 25 if possible. A further suggestion 
was that when this value was exceeded special atten- 
tion should be given to the factors that insure lateral 
rigidity. (Reference 10, p. 16, and 1923 annual 
report, p. 390.) As a result of the present experiments, 
the Forest Products Laboratory has learned what 
factors are involved in the lateral buckUng load and has 
concluded that no arbitrary ratio for the moments of 
inertia can properly be set and that such a method of 
design should not be used. 

In previous tests it was practically impossible to 
prevent the buckling of I beams having a moment-of- 
inertia ratio of 39. In the panel with the 1 by 6 inch 
beams just mentioned, for which the moment-of- 
inertia ratio is 36, the maximum moment w^as approxi- 
mately 89 per cent of the moment that would have 
been required to cause failure had bending been 
confined to a vertical plane, and even this percentage 
value could not have been obtained if it had not been 
for the excess strength of the ribs. In the third test of 
the last panel, which had K by 4 inch beams and for 
which the moment ratio is 64, the maximum load was 
approximately 76 per cent of the load required to 
cause failure had bending been confined to a vertical 
plane. 

In all of the recent tests it is probable that the 
beams were receiving less lateral support than the 
beams in an ordinary wing panel would receive and the 
end fixity was less than that which obtains in the usual 
drag bay. With a more or less rigid nose, such as one 
of plywood or metal, and ribs slightly over strength, 
beams with moment-of-inertia ratios considerably in 
excess of 25 can be counted upon for their full bending 
strength. 

CONCLUSIONS FOR PART II 

Deep beams may fail through buckling laterally and 
twisting at loads much less than those calculated by 
means of the usual beam formula. 

There is for each fixity and loading condition a 
critical lateral buckling load for a deep beam just as 
there is a critical load for a column. 

A mathematical analysis of the problem for various 
loading and fixity conditions leads to formulas that 
contain the dimensions of the beam, the modulus of 
elasticity along the grain, the modulus of rigidity in 
torsion, the span, and a constant depending upon the 
loading and fixity conditions. 

Experimental results confirm the practical appU- 
cabihty of these formulas. 

When one spar of an airplane wing or other panel 
is subjected to an axial load and the other spar and the 
ribs are not loaded, the lateral rigidity of the whole 
combination is the sum of the lateral rigidities of the 
two spars. 
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When two deep beams fastened together with ribs 
are subjected to bending, lateral buckling of the beams 
may or may not be prevented. When one or both of 
such beams are heavily stressed and in need of lateral 
support, the ribs, if they are not stronger than is neces- 
sary to carry the load upon them, can not carry the 
extra load that is induced by the tendency of the 
beams to buckle. 

A fairly rigid nose and ribs slightly overstrength will 
permit the use of aircraft wing beams that have a rela- 
tively large ratio of moment of inertia about the prin- 
cipal horizontal axis to that about the principal vertical 
axis. 

No arbitrary moment-of-inertia ratio can be used 
with certainty. Each particular case must be studied 
individually and lateral support must be provided in 
accordance with the tendency of the beam to buckle 
laterally rather than to bend in a vertical plane. 

This investigation was undertaken as a study in air- 
craft design. The conclusions, however, are of general 
application, even though some of them for convenience 
are worded as if they applied only to aircraft. 
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APPENDIX 



A MATHEMATICAL STUDY OF THE ELASTIC STABILITY OF THIN, OUTSTANDING FLANGES 

UNDER COMPRESSION 



INTRODUCTION 

III discussing the stability of a column or otlnr 
compression member having one or more thin, out- 
standing flanges, it is necessary to consider not only 
^ the conditions for the stabiUty of the 

column as a whole but also the stability 
of the flanges themselves. The problem 
of the stability of such a flange is essen- 
tially that of the stability of a rectan- 
gular plate simply supported along the 
ends to which the load is applied, frre 
along one of the other edges, and on the 
remaining edge either simply supported, 
imperfectly fixed, or perfectly fixed, de- 
pending upon the nature of the section. 
Timosheno has discussed this problem 
in considerable detail for plates of 
isotropic material. (References 17 and 
21.) In the following appendix his 
methods will be extended to plates com- 
posed of a nonisotropic material, such 
as wood, which will be considered 
to have three mutually perpendicular 
planes of elastic symmetry. His anal- 
ysis for isotropic plates will also be 
summarized and some further conclu- 
sions drawn. 



Figure 18.— A rec- 
t a D g u 1 a r plate 
under a uniform 
compressive load 
on two opposite 



EXACT METHOD; BASE OF FLANGE PERFECTLY FIXED 

DIFFERENTIAL EQUATION FOR THE DEFLECTION OF A FLANGE 
OF NONISOTROPIC MATERIAL UNDER A COMPRESSIVE LOAD 

A plate of thickness A, Figure 18, is considered to 
he in the XF-plane and to be bounded by the hues 
j = 0, x = tt, y = 0, and y = h. Uniform compressive 
loads P per unit length of edge, parallel to the ^-axis^ 
are applied to the edges a- = 0 and x' = a, which are 



simply supported. The edge y = b is free while the 
edge 2/ = 0 is either simply supported, partially fixed, 
or perfectly fixed. 

The case in which the edge y = 0 is perfectly fixed, 
a case which rarely or never occurs in practice, is first 
treated for both isotropic and nonisotropic material, 
making use of the differential equation for the deflec- 
tion of the plate from its plane and of appropriate 
boundary conditions. A simpler approximate method 
based on energy considerations is then applied to the 
same case and the results are compared and found to 
check in a satisfactory manner. The approximate 
method is then applied to the case in which the edge 
in question is only partially fixed, the case in which 
the edge is simply supported appearing as a h mi ting 
form of partial fixity. 

The differential equation satisfied by the deflection 
w is obtained from the following differential equations 
connecting the stress resultants T, Sj and A^^ and the 
stress couples G and H acting upon an elementary 
portion of the plate with edges dx and dy. (Reference 
8; art. 326, equations (24), (25), (26), and art. 331, 
equations (45) and (46).) The notation used is that 
of Love. (Reference 8, art. 294.) 

dT\_ dSo dhv_ d'w 
dx dy' ^'dx' 
d-iv 



dx ^ dy dxdy 



dN^.dN, 
dx dy 



dx' 



^ dxdy 
d'w 
dy 
d^w 



(15) 



d^w ^ 
'^'dxdy^^' dxdy 



dhv 

dy' 



4-Z' = 0. 



d^w 
dydx 



-Go 



dH, 
dx 
dG, 
dx 

dho_ 
dxdy 



+ 



dG, 
dy 
dH2 
dy 



(16) 



d'w rj 



dhv 



+ Si + S, = 0. 





Figure 19— (A) Stress resultants and (B) stress couples on an element of a plate 
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In equations (15) and (16) X' =^Y' = Z' =M' = N' = 0, 
since the components of the external force per unit 
area and of the external couple are zero. 

To calculate Ti it is necessary to 

express the components of stress Xx XyXa 

terms of the deflection w and the elastic constants- 
(Figure 20.) The displacements u and v are given 
with sufficient accuracy by 

dx 
dw 



u = 



(17) 



V= —2 



dy 



The components of strain are 



— z 



du 
dx 

_dv _ _ 



_du dv_ _ d^w 
^''~dy^dx~ ^^dxdy' 



dx' 



(18) 



For a more extensive discussion of the components of 
strain, see article 329 of reference 8. 

Assume that the material of the plate, wood, has 
three mutually perpendicular planes of elastic sym- 
metry. (Reference 8, arts. 110 and 111.) Denote by 

0 




Figure 20.— Components of displacement in terms of 
deflection 

Ey, and Ez Young's moduli in the directions x, y, 
and z, respectively, by a^y Poisson's ratio associated 
with contraction parallel to the F-axis and stress 
parallel to the X-axis, and by Mzv the modulus of rig- 
idity corresponding to the directions x and y. Tne 
stress components Xxy Yy, and Xy are then given by 



X, 



Ex 



1 ^xy^ yx 



{Cxx "i" O'yx^yy^ 



(19) 



1 



Xy fXxy ^xy 2fixyZ 



dxdy' 



By definition 



Jh/2 
Xxzdz, 
-hl2 



(Reference 8, art. 294.) Then 



^\dx' 



+ (7„ 



'dy')' 



(20) 



where 



In like manner 



\2{\— ajyOyl) 

d'w 



~ ~^\'dy^'^'''''~d?)' 



where 



12(1 -<Tiy<Tyx) 

Further, from their definitions, 



where 



M=Hxy-^ 



dxdy 



(21) 
(22) 

(23) 

(24) 
(25) 



In the last of equations (16) the quantities G^i, 6^2, 
H^, and H2, which are expressed by (20), (22), and (24) 
in terms of second partial derivatives of Wy are each 
multiplied by second derivatives of w. Each of these 
derivatives may be considered small and the product of 
two of them negligible. It follows that 



Si — — S2 



(26) 



From the first two equations (16) and equations 
(20), (22), and (24) it is found that 



AT /^^^ . \ TIT ^'^ 



iV2 = 



■ ^2 ^Qy^ + cTxy Q^2Qyj- 



(27) 



M 



It is clear from their definitions and (19) that Si and 
^2 are small. (Reference 8, art. 294.) Also from its 
definition and equation (18) To is small. Equation 
(26) and the first two of equations (15) are satisfied 
approximately by taking 



and 



Si=S2=T2 = 0 
Ti = constant = — P 



(28) 
(29) 



where P is the load per unit length of the loadeded ges. 
The third of equations (15), on making use of (27), 
(28), and (29), then gives the differential equation of 
the plate : 



Or 



d*w 



d*iv^p d'w 



where 



K= 



dx'dy'~^^'dy'^' dx' 
Di<7yi + D2ajy + 2M 



= 0 



(30) 



(31) 
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BOUNDARY CONDITIONS 

On the simply supported edges x = 0 and x = a 
Figure 18, 

w = 0 (32) 

and 

Gi = 0. (32a) 



The last condition requh-es that 



dx 



,2 + O'yx n-.2 



On the fixed edge, y = Oj 

w = 0 



and 



dw 

dy' 



0. 



(33) 

(34) 
(35) 



On the free edge, y = by 

G2 = 0 



and 



dH2 



dx ~ vat/ " dx'dy) dx'^dy 
Rewriting these conditions for the edge y = h, 



2M 



-=0. 



dy 



72 + 



= 0 



and 



ivhere 



that is, 



dy3 + (^ dx'dy' 



0, 



(36) 
(37) 



2 — <r = crxu + - 



2M 



(2— aiy)Ey—4lX:ru (1 ~ ffjrKffyT) 



(38) 



SOLUTION OF THE DIFFERENTIAL EQUATION 

Conditions (32) and (33) are satisfied by 



It; = sin ??^y (y) = sill \xf(y). 



(39) 



It will be convenient to replace rmr/a by tt/c, for if the 
flange breaks up into more than a single half wave 
each portion of length a/m = c may be considered as a 
plate of length c simply supported at its ends. We 
shall accordingly interpret X as given by the equation 



X = - 



where c may be either the entire length of the flange or 
a portion of this length, as circumstances require. 

In accordance with (SO) f(y) in (39) must satisfy an 
ordinary linear difl'erential equation of the fourth order. 
Its solution can be written 

fiy) = C7ie-«^ + C2C«^ + cos Py + d sin py, (40) 



where 



(41) 



Conditions (34) and (35) are satisfied if the constants 
in (40) are so related that we may write 

f{y)=A (cos I3y~ cosh ay)^B (sin /37/-^sinh ay). (42) 

The substitution of (39) combined with (42) in the 
conditions (36) and (37) leads to the equations: 

A[{0'^<T,yK') COS fihA- {a^-(j,y\^) cosh + 

5[(i32+(7,,X-) sin ^6 + ^(a--o-,,X2) gj^h ah] = 0 (43) 

and 

.4[iS(/5- + 2X2-(rX2) sin fib - a{a^ -2\^ -\- a\^) sinh a6] + 
J5[-i3(iS- + 2X2-(7X2) cos i36-/3(a— 2X2 + (tX2) cosh ah] = 0. 

(44) 

In (44) note that after some reduction 

/32+(2-(7)X2 = a2-(r^,X^ 
a^-{2-(T)\^ = 0'^cT,y\\ 

In this reduction the following relations were used: 

2M . 

2 — (J= (Txv^ *Ti~ ' and Ey(Txy = Ex(Tyx' 



(45) 



(Reference 1, p. 104.) Using the abbreviations 

the equations (43) and (44) can be written in the form 

A[t cos 13b + s cosh ab] + B[t sin ^b-^~s sinh ab] = 0 

(46) 

A[ps sin 13b - sinh a6] + B[- 13s cos i86 - cosh aft] =^ 0. 

In order that solutions of the system (46) other than 
.4 = 0 and B = 0 may exist, that is, that a solution 
different from zero of the differential equation (30) of 
the form (39) may exist, it is necessary and sufficient 
that the determinant of the coeflBcients of A and B in 
(46) vanish. The result of equating the determinant 
to zero is, after some reduction, 

2^6'+ (f + s'') cos ^b cosh ab = l — ' J sin fib sinh ab. 

(47) 

Multiplying this equation through by b\ the terms can 
be arranged so that a and fi occur only in the combina- 
tions ab and fib. We then write (equations (41)) 

ab=uuv^-vy 

(48) 

fib = UUV-V)\ 



where 
and 



V=^^ (49) 
U= j^^j^ W - ^. A.) + D,PW] . (50) 
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GENERAL EXPRESSION FOR CRITICAL STRESS 

By assigning a value to the ratio cjh the quantity 
V is determined. The corresponding value of U can 
then be found by solving equation (47). The value 
of the critical stress 

P 

^ = 1 

corresponding to this value of c/6 can then be found 
from equation (50). From (50) it follows that 

_r K J, 1 /K' fi' 

- \_EJ^ ^'"^ c'EJ^' V A " Vj^ P ■ 



Or 



(51) 



ELASTIC CONSTANTS OF SPRUCE 



The elastic constants to be used in the computa- 
tion depend upon the orientation of the planes of 
elastic symmetry of the wood in the plate. It will 
be assumed throughout the discussion that the grain 
of the wood is ])ai-allel to the .Y-axis, the direction 
in which the compressive load is applied. Two cases 
for the direction of the growth rings of the wood will 
be considered, one in which the rings are perpendicu- 
lar to the faces of the plate and another in which 
they make an angle of 45^ with the faces. 

In the first case (fig. 21) Young's moduli E^y Ey, and 
^ Ez are equal to El, Er, and Et, 

respectively, the subscripts Z, R, 
and T denoting the longitudinal, 
radial, and tangential moduli, re- 
spectively. The values for these 
and other elastic constants for 
spruce were taken from a report of 
the British Aeronautical Research 
Committee. (Reference 1, p. 105.) 
The values are: 



Figure 21.— The cross 
section of a quarter- 
sawn flange 



Er. 


-1. 95X10^ 


o'lr = 


0. 45 


Er 


= 0. 13X10^ 


(^LT = 


0. 539 


Et 


-0. 07X10^ 


(TrT — 


0. 559 ■ 




-0. 104 XIO^ 


(^RL = 


0. 03 


f^LT 


= 0. 072 XIO'^' 


(^TL = 


0. 0194 




= 0. 005X10^ 


o'tr'^ 


0. 301 



In the second case, when the growth rings make an 
angle of 45° with the faces of the plate (fig. 22), the 



elastic constants Ey /jLjry can })e com- 
puted from those just given by the following forniuhis: 



Ey 4En4Er'^4nRr 2Er 

Ey 

^l^LTl^LR 



(Reference 8, art. 111.) It is 
then found that 



E, 



^0.01875X 10^ 
=0. 00475 




a,y = 0. 494 



tx^y = 0. 0851 XIO^ 



Figure 22.— The cro.ss section 
of a wood flange the growth 
rings of which make an angle 
of 45° with the faces 



CRITICAL STRESS FOR A FLANGE OF SPRUCE 

Yalues of k in equation (51), the equation for critical 
stress, which result from solving equation (47) for the 
cases of growth rings perpendicular to the faces of the 
flange and at 45"" to the faces, are given in Tables YI 
and YII, respectively. 

TABLE VI 

THEORETICAL CONSTANTS FOR FLANGED COM- 
PRESSION MEMBERS OF SPRUCE HAVING THE 
GROWTH RINGS PERPENDICULAR TO THE FACES 
OF THE FLANGE, CALCULATED BY THE EXACT 
MATHEMATICAL METHOD 



c/b 


U 


k 


3. 40 


15. 42 


0.228540 


3.30 


15. 13 


.228345 


3.25 


14. 95 


.228057 


3.20 


14. 82 


.228376 


3. 10 


14.55 


.229166 


3.00 


14.30 


.230677 



The minimum critical stress for growth rings perpen- 
dicular to the faces of the flange occurs when the half 
wave length is 3.25 times the outstanding width of the 
flange. This critical stress is equal to 0.228 EJiyb\ 
Ordinarily the length of the colmnn is such that the 
flange can not break up into segments the length of 
which is exactly 3.25 times the outstanding width. 
Under such a condition the stress will be increased as 
the values in the table indicate. Considerable increases 
would be found for considerable departures from the 
optimum value of the ratio c/b. Such departures occur 
only when the column is so short that its length is less 
than two or three times the optimum half wave length. 
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TABLE VII 

THEORETICAL CONSTANTS FOR FLANGED COM- 
PRESSION MEMBERS OF SPRUCE HAVING THE 
GROWTH RINGS AT AN ANGLE OF 45° WITH THE 
FACES OF THE FLANGE, CALCULATED BY THE 
EXACT MATHEMATICAL METHOD 



c/b 


U 


k 


5. 34 


14.83 


0. 117308 


5.20 


14. 78 


.117119 


5. 10 


14. 76 


. 117113 


5.00 


14. 74 


.117118 


4.90 


14.73 


. 117227 


4.80 


14. 75 


. 117539 



Consideration of Tables VI and VII shows that the 
theoretical critical stress is considerably less when the 
growth rings make an angle of 45° with the faces 
of the flange (fig. 22) than when they make an angle 
of 90° (fig. 21). The chief factor in determining the 
variation in the critical stress with variation in the 
angle between rings and faces is the ratio EyjE^. 
This ratio is nearly constant when the angle made by 
the rings with the faces of the flange lies between 20° 
and 70°, and hence the results for rings at an angle 
of 45° may be taken to apply over this range. When 
the rings are parallel to the faces of the flange, 
however, the minimum critical stress is found by an 
approximate method given later in this report to be 

- 0.164 E.h'lh' 
for a flange with a perfectly fixed edge. This critical 
stress is intermediate between those for flanges with 
the rings at angles of 45° and of 90° with the faces. 

The theoretical critical stress for a flange with a 
perfectly fixed edge is not attained in practice because 
the condition of perfect fixity at the base of the flange 
is not realized. Later in this report it will be pointed 
out more in detail that as the fixity at the base of the 
flange decreases the variation of the critical stress with 
inclination of growth rings becomes smaller and 
ultimately, as the fixity continues to diminish, the 
critical stress for a flange with rings parallel to its 
faces becomes less than that for a simihir flange with 
rings at 45°, which in turn is always less than that for 
a similar flange with range at 90°. 

DIFFERENTIAL EQUATION AND BOUNDARY CONDITIONS FOR A 
FLANGE OF ISOTROPIC MATERIAL 

The preceding analysis is an extension to flanges 
of nonisotropic material of the method that Timo- 
shenko used in discussing flanges of isotropic material. 
(Reference 17, p. 350.) When the material is isotropic 
the differential equation (30) becomes 



where 



(7=: 



Eh? 



(52) 



(53) 



■l2(l-<r^) 

The boundary conditions are given by equations (32) 
to (37) after a^y and (r„ have been replaced with a. The 



differential equation and the boundary conditions are 
then those used by Timoshcnko. The critical load is 
determined by solving equation (47) where i and s 
are given by (45) with (rx^ = o- and where a and /3 are 
given by (48) with 



and 



(54) 



(55) 



CRITICAL STRESS FOR A FLANGE OF ISOTROPIC MATERIAL 

The values of U corresponding to various vakies of 
the ratio cjb as calculated by Timoshenko are given in 
Table VIII for flanges of isotropic material. In the 
third column of this table appear the values of k in 
the formula 



V 



where p is the critical stress. This formula is obtained 
at once from equation (55) by noting that 



^^h ^=12(1 



EK^ 



V) 



In the computations cr was taken as 0.25. 

TABLE VIII 

THEORETICAL CONSTANTS FOR FLANGED COM- 
PRESSION MEMBERS OF ISOTROPIC MATERIAL, 
CALCULATED BY THE EXACT MATHEMATICAL 
METHOD AND WITH POISSON'S RATIO TAKEN 
AS 0.25. 



cjb 


U 


k 


1.0 


16. 76 


1.490 


1. 1 


15.41 


1. 370 


1.2 


14.47 


1.286 


1.3 


13.88 


1. 234 


1.4 


13. 45 


1. 196 


1.5 


13.20 


1. 173 


1.6 


13. 13 


1. 167 


1.635 


13. 1 1 


1. 165 


1. 7 


13. 15 


1. 169 


1.8 


13. 24 


1. 177 


1.9 


1.3. 43 


1. 194 


2.0 


13. 67 


1.215 


2.2 


14. 35 


1.276 


2.4 


15. 21 


1. 352 



In Table VIII the critical stress is least when the 
half wave length is equal to 1.635 times the width of 
the outstanding flange. If a, the total length, is either 
less than 1.6356 or somewhat greater than this amount 
the critical stress will be greater, as Table VIII shows. 
As a increases toward twice the ideal half wave length 
the critical stress begins to diminish, reaching the 
same minimum value at a = 3.276 as at a = 1.6356. 
When the column is long in comparison with the 
width of the outstanding flange (the length three or 
more times the width) the flange will break up into 
waves the half length of which is approximately 
1.6356, and the critical stress will then differ but little 
from that for this ideal half wave length. 
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APPROXIMATE METHOD 

DISCUSSION 

Approximate results were obtained by Timoshenko 
with a method that is based upon energy relationships 
and that is an important extension of a method used 
by Bryan. (References 3, 4, 5, 19, 20, and 21.) The 
deflection of the plate (fig. 18) is expressed as a sum of 
terms of the form 

w = Ai(l>i(:x, y)-\rA2<t>2(x, y) + . . . . (56) 

the functions </>2, being chosen to satisfy 

the boundary conditions as nearly as possible and the 

coefficients yli, being arbitrary. This 

expression for the deflection w is then substituted in 
the integral representing the energy of deformation of 
the plate. The result is a function of the arl)itrary 

constants Ai, An, The energy is then 

equated to the work done by the compressive load P 
per unit length acting on the edges x = 0 and x = a» 
The result is an equation that can be solved for P 

Ml terms of the arbitrary constants Ai, Aoy 

The ratios Ao/A^ Az/Aij are then chosen in 

such a way as to make P a minimum. If the resulting 
stress, 

P 

where k is the thickness of the plate, is less than the 
stress for primary failure of the column of which the 
plate is a member, the plate will fail by buckling at the 
critical stress p. For a full discussion of the method, 
with examples of its application to simple cases, see 
Timoshenko 's paper. (Reference 19.) 

The energy of deformation of the plate, under the 
assumption that the stress components Xzy Yzj and 
Zz are negligible, is given by 

h_ 

Substituting the values of the strain components given 
in (18) and those of the stress components given in 
(19) for nonisotropic material having three mutually 
perpendicular planes of elastic symmetiy, the result is 

^^lijo Jo |r-<r..<r„I^'^(^x^) +^'"(5^) 

+ 2E,<tJ^, g5] + 4Mx.(at4)V2/ dx. (58) 
For isotropic material this boconies 



The work done by a compressive load P per unit 
length of edge, applied to the edges x = 0 and x = a 
(fig. 18) is given by 

In what follows, the integrations with respect to x 
in (58), (59), and (60) will be performed between the 
limits 0 and c, where c is the half wave length of the 
deformed surface. In certain cases c will be equal to 
a, while in others it will be a fractional part of a. 

BASE OF FLANGE PERFECTLY FIXED 

The assumed deflection (equation (56)) will be 
taken as 

w={Ax {QbY - 462/' + y') + A2{y^ - 106 V 

+ 206y)} siny. ^^^^ 

The functions of 7/ in the first and second terms of 
(61) represent respectively the deflection of a canti- 
lever fixed at the end y = 0 under a uniform load and 
under a load that is proportional to y. Timoshenko 
in treating the isotropic plate by this method chose 
other functions. (Reference 21, p. 405.) It is not 
apparent that either choice possesses any particular 
advantages over the other. 

Flange of nonisotropic material. 

Entering (61) in (58) it follows that for nonisotropic 
material 

in which, letting 

P = p» (63) 

(/o = 2.31l7r^4-7r2 
[41.15 ^-y^^-£^y^y^^ _ 3 .432(r J p + 28 .8 

(/, = 16.7887r^+7r2 

[303 .4 fL^ySk^^y^l _ 26.30cr,,] p + 208.3 |j p' (64) 
(/2 = 30.4887r*+7r- 

[559.7 ^'^y^^-^ yS-y^ - 50. 16(7,.]p + 377.2 1^ 

and 

2 = ^6. (65) 
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From (60) and (61) it follows that 



T=f^A,'h\c,-\-c,z^C2z') 



(66) 



where 



Co 



= 2.311, ci = 16.788, and C2 = 30.488. (67) 



Equating T and V as given by (62) and (66) and solving 
forp = P/A, 



I27r2(l - (7:,y(T,^)p\Co-^rC,Z + CoZ'') ""b^' 



(68) 



The critical stress p will be a minimum if z is the 
larger of the roots of 

dz 



Equation (68) may be written 



(69) 



The calculation outlined assumes the ratio c/b to be 
given and determines the critical stress for this ratio. 
By calculating k for a series of ratios c/b the ideal half 
wave length is found as that which makes the critical 
stress a minimum. 

In Tables IX, X, and XI, the values of k for cer- 
tain values of the ratio c/b are given for flanges of 
spruce, the growth rings being respectively perpen- 
dicular to the faces of the flange, inclined to them at 
an angle of 45°, and parallel to them. The elastic 
constants for spruce given earlier in this appendix 
were used in the calculations. For rings parallel to 
the faces, we note that 

Ej,= l. 95X10^ 

Ey = 0. 07X10^ 

fx,y = 0. 072X10^ 

(Txy = 0, 539 

ayx = 0. 0194. 
(Reference 1, p. 105.) 

TABLE IX 

THEORETICAL CONSTANTS FOR FLANGES OF 
SPRUCE, UNDER LONGITUDINAL COMPRESSION, 
THAT HAVE THE GROWTH RINGS PERPENDICU- 
LAR TO THE FACES OF THE FLANGE, CALCULATED 
BY THE APPROXIMATE METHOD 



c/6 


k 


3.2 
3.3 
3.4 


0. 228356 
. 228256 
.228719 



TABLE X 

THEORETICAL CONSTANTS FOR FLANGES OF 
SPRUCE, UNDER LONGITUDINAL COMPRESSION, 
THAT HAVE THE GROWTH RINGS AT AN ANGLE 
OF 45° WITH THE FACES OF THE FLANGE, CAL- 
CULATED BY THE APPROXIMATE METHOD 



c/b 


k 


5. 1 
5.2 

5.3 


0. \ \m2 

. 1 16806 
. 117021 



TABLE XI 

THEORETICAL CONSTANTS FOR FLANGES OF 
SPRUCE, UNDER LONGITUDINAL COMPRESSION, 
THAT HAVE THE GROWTH RINGS PARALLEL TO 
THE FACES OF THE FLANGE, CALCULATED BY 
THE APPROXIMATE METHOD 



cjb 


k 


3.7 
3.8 
3.9 


0. 16399 
. 16381 
. 16389 



The results agree remarkably well with those given 
in Tables VI and VII as the result of more exact 
analysis. 

Flange of isotropic material. 

After substituting the assumed deflection (61) in 
the integral (59) for the energy of deformation of the 
flange in the case of isotropic material and equating 
T and V as given by (59) and (60) it is found that 

^ ^ 12^\l-a^)p\7o + c,z + C2z'rb^ ^ ' 

where 

P 

P=h 

and expressions for doy rfi, (1-2 are found from (64) by 
writing 

The quantities Co, Ci, and Co have the values given by 
(67). 

If equation (70) is written in the form 

V = kE]i, (71) 



the value of the minimum k for a given value of the 
ratio c/b can be calculated as with nonisotropic mate- 
rial. A few values in the vicinity of the half wave 
length for which the critical stress is a minimum are 
given in Table XII; Poisson's ratio o- was taken as 
0.25. 
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TABLE XII 

THEORETICAL CONSTANTS FOR FLANGES OF ISO- 
TROPIC MATERIAL UNDER LONGITUDINAL COM- 
PRESSION, CALCULATED BY THE APPROXIMATE 
METHOD AND WITH POISSON'S RATIO TAKEN AS 
0.25 



c/6 


k 


1. 65 
1.66 
1.67 


1. 16434 
1. 16390 
1. 16407 



The minimum values of k in Table XII differ from 
those of Table VIII by a small fraction of 1 per cent. 
The half wave lengths at which the minimum critical 
stress occurs differ by about 1.5 per cent. Plotting 
the curve connecting critical stress and half wave 
length in the vicinity of the minimum critical stress 
will show that this difference has little significance. 
For steel, with Poisson's ratio taken as 0.3, a similar 
calculation gives a minimum k of 1.1592 correspond- 
ing t(5 a value of c/b of 1 .60. 

BASE OF FLANGE IMPERFECTLY FIXED 

Discussion. 

The condition of perfect fixity assumed in the pre- 
ceding sections of this report for the edge of the flange 
y = 0 (fig. 18) is probably never realized. This is due 
to two circumstances, which will be considered sepa- 
rately. Both result from the moment induced at the 
edge 2/ = 0 by the deformation of the outstanding flange 
bounded by this edge. This moment causes twisting 
of the whole cross section of the column and it also 
causes elastic giving of the material along the junction 
of the base of the flange and the body of the column. 
Both of these phenomena, twisting of the section and 
elastic giving at the base of the flange, are accompanied 
by a change in the inclination of the flange at its base 
from the value zero required by the condition of per- 
fect fixity. The twisting phenomenon is easily ex- 
pressed in terms of the torsional rigidity of the section. 
The elastic giving appears to involve factors that are 
best determined experimentally. 

Effect of twisting of column. 

We proceed to calculate the effect of the twisting of 
the column induced by the moments acting along the 
edge y = 0. (See Timoshenko. Reference 21, p. 400.) 

Let </) denote the angle of rotation of a cross section 
the abscissa of which is x. If elastic giving of the 
material is neglected for the present, 



The torsional couple in any section is then 



(72) 



dx 



^^i^dxdy)y„o 



where G is the modulus of rigidity of the material and 
Kis the torsion constant of the section. (Reference 22, 



p. 11, and 1929 annual report, p. 681.) 
appHed per unit length is then 



The couple 



ydx^'dyjy^o 



In applying the approximate method, the strain 
energy resulting from the twisting of the column (in 
whole or in segments) should be added to the strain 
energy of deformation of the corresponding portion of 
the outstanding flange. The strain energy per half 
wave length c, resulting from twisting, is 



If 



'i^=J(y) sin 



TTX 



V,- 



To apply the approximate method let 
w = ^Ay + ^1 (l - cos -^^'j sill 



irx 



(73) 



(74) 



If ^1 = 0 the edge y = 0 is simply supported. If ^1 = 0 
the edge y = 0 is fixed. Hence, by allowing AJA to 
vary from zero to infinity, all conditions on the edge 
y = 0 intermediate between those for an edge simply 
supported and those for one perfectly fixed can be 
satisfied by a deflection in the form given by (74). 
It follows from (73) and (74) that 



A'c 



where 



4 



(75) 
(76) 



In calculating Vi for a column of spruce the modulus 
of rigidity G may be taken as the mean modulus that 
would be given by a torsion test on a cylinder of circular 
section. This value may be conveniently taken as 
Young's modulus in the longitudinal direction divided 
by 15.6. (Reference 22, pp. 21 and 24, and 1929 
annual report, pp. 091 and 094.) 

Flange of nonisotropic material. 

For nonisotropic material, such as wood, with three 
mutually perpendicular planes of elastic symmetry it 
follows from (58) and (74) that 

E^h'cT' f . , r_6^ . 4Mx.(l-cr,,(r,,) J_-| 



48(1 — (Txydyx) 



+ 



irc^ \ IT/ 



+ 



2c2, 
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From (60) and (74) 



From 
using 



A, , P 



(77) 



it follows that 



or 



(78a) 



where for convenience the following notation has been 
used: 



, .,GK b 
= 12(1 - o-j„(T„r)ir--^^2/,> 



(79) 



, _ 1 . 4;Uxy(l - q-iyO-i,i) p 



rf. = l-- + ^-^<l-^)p+^^^"^T^ V (80) 



(/2 = 



3 4 

9. w 



Er 



2"^ £'.32 



1 ,4,8 , 34 

Co=o' Ci = l + -2' '^"dC2=o""~" 



(81) 



Flange of isotropic material. 

For isotropic material equation (78) with appropri- 
ate values of do, du ^^2, Co, Ci, and C2 becomes that given 
by Timoshenko, to whom the choice of the form (74) 
for the deflection is due. (Reference 21, p. 401.) 

The vahies of the coefficients d^, di, and do are 

; 1 , 2(l-(r) ep' 

, 3 4 , /I crV 
^*=2~^ + ''U~'^/ 32 



where 



ion 2^ .GK h 
e=12(l-<r)7r -g 



(83) 



The constants cn, Cx, and C2 are unchanged. 



The critical stress p is given by 



12(1 -<r^)pV 



or 



(84) 
(84a) 



Application of formulas. 

E(l nations (78) and (84) are of the same form as (68) 
and should be used in the same way. For a given e 
and a series of values of the ratio c/6 a series of critical 
stresses p are determined corresponding to a suitable 
value z. The ratio c/b associated with the minimum 
critical stress (if there is a minimum) determines the 
half wave length c that is ideal for the value of e under 
consideration. 

For the study of a given column it is more conven- 
ient to proceed in another way. The first step is to 
construct a table giving k in the formula for the 
critical stress as a function of the fixity coefficient e, 
for each of a series of values of the ratio c/b of the half 
wave length to the width of the outstanding flange. 
Table XIII was constructed in this way for flanges of 
si)ruce and Table XIV for flanges of isotropic material. 
The results in these tables are also shown in the curves 
of Figures 23 and 24. 

The use of these curves in studying a given column 
is discussed in a later section of this appendix. In 
interpreting the curves, it must be borne in mind that 
the fixity coefficient e depends upon the half wave 
length c and the outstanding width b as well as upon 
the thickness h and the torsion constant K. 

TABLE XIII 

THE COEFFICIENT k IN EQUATION (78a) FOR A 
FLANGE OF SPRUCE HAVING GROWTH RINGS AT 
AN ANGLE OF 45° WITH THE FACES 



c/b 




k 


c/b 


€ 


k 




0. 10 


0. 868 


3 


0. 10 


0. 154 


1 


.05 


.866 


3 


.05 


. 147 


1 


.03 


.865 


3 


.03 


. 142 


1 


.01 


.864 


3 


.01 


. 136 




.00 


.863 


3 


.00 


. 132 


5 


.10 


.113 


7 


.10 


.117 


5 


.05 


. 102 


7 


.05 


. 102 


5 


.03 


.094 


7 


.03 


.091 


5 


.01 


.083 


7 


.01 


.073 


5 


.00 


.075 


7 


.00 


.059 


9 


.10 


.136 


12 


.10 


. 181 


9 


.05 


. 115 


12 


.05 


. 149 


9 


.03 


.100 


12 


.03 


. 124 


9 


.01 


.074 


12 


.01 


.079 


9 


.00 


.053 


12 


.00 


.049 


15 


.10 


.241 


20 


. 10 


.375 


15 


.05 


.195 


20 


.05 


.298 


15 


.03 


.159 


20 


.03 


.237 


15 


.01 


.098 


20 


.01 


.133 


15 


.00 


.047 


20 


.00 


.046 



1 
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TABLE XIV 

THE COEFf^ICIENT k IN EQirATION (84a) FOR A 
FLANGE OF I80TR0PIC MATERIAL THAT HAS 
A POISSON'S RATIO OF 0.25 



C/0 


c 


f. 

K 


1 


0.50 


1. 287 


1 


. 10 


1. 273 


1 


.00 


1.269 








3 


.50 


.616 


3 


. 10 


.522 


3 


.05 


.509 


3 


.01 


.498 


3 


.00 


.496 


5 


.50 


.750 


5 


. 10 


.503 


5 


.05 


.469 


5 


.02 


.448 


.5 


.00 


.434 


10 


.50 


1.640 


10 


. 10 


.676 


10 


.05 


.543 


10 


.00 


.409 



c/6 



0.50 
. 10 
.05 
.01 
.00 

.50 
. 10 
.05 
.03 
.01 
.00 

. 50 
. 10 
.05 
.03 
.01 
.00 



0. 673 
.628 
.622 
.617 
.616 

.660 
.498 
.476 
.467 
.458 
.4.54 

1.026 
. .549 
.484 
.4.58 
.431 
.418 



Flange with a simply supported edge, the limiting case 
as the fixity coefficient approaches zero. 
As the fixity coefRcient e ap])roaches zero in equations 
(78) and (84) it is found that the value of z correspond- 
ing to a minimum value of y approaches zero. This 
should be so for as e approaches zero the edge y = 0 
becomes more and more nearly simply supported. 
The ratio of Ai to A in (74) will then approach zero. 
By equation (77) this implies that z approaches zero, 
as just noted. 

Accordingly the limiting critical stress as e ap- 
proaches zero is found to be 



12(1- 



1 _^ fJ-xy 



(85) 



by setting € = 0 and £: = 0 in (78) and (SO). The vahies 
of k given by this formula for a simply supported edge 
agree well with those of Table Xlll for the fixity 
coefRcient e = 0. As p becomes large p decreases to the 
limiting value, 



Using the elastic constants for spruce having the 
growth rings at an angle of 45° with the faces of the 
flange, (86) becomes 



P = 0.04:4: E, 



(87) 



If the growth rings are perpendicular to the faces of the 
flange 

2? = 0.053 eJ^^. (88) 
while if they are parallel 

p = 0.037 (89) 



Thus for a flange with a simply supported edge the 
critical stress is less when the growth rings are parallel 
to the faces of the flange than when the rings are 
inclined to them at an angle of 45°. For a flange with 
a perfectly fixed edge, on the contrary, the critical 
stress was found to be less when the rings are inclined 
to the faces at an angle of 45° than when they are 
parallel to them. The relative variation of the critical 
stress with inclination of the rings is less for flanges 
with simpl}^ supported edges than for those with 
perfectly fixed edges. 

In practice, the fixity at the bases of the flanges is 
small. Consequently the variation of the critical 
stress with the inclination of the growth rings may be 
expected to be similar to that for flanges with simply 
supported edges. 

From this point on the discussion will be limited to 
flanges with growth rings at an angle of 45° to the 
faces. The results may be considered to be applicable 
to flanges with rings at any inclination except for the 
extreme cases of rings nearly parallel to the faces or 
nearly perpendicular to them. In the first case the 
calculated critical stress should be reduced somewhat, 
while in the second it should be increased somewhat. 
These formulas hold for long flanges. For short ones 
the effect of the first term of (85) must be included. 

For isotropic material the equations corresponding 
to (85) and (86) are 



V 



and 



-[ 



12(1 -<7^)p^ 2(1 



1 



With 



2(l + <r)""62 
= 0.25 equation (91) becomes 



■p = QAE 



and with <r = 0.3 



2> = 0.385£'p- 



(90) 
(91) 

(92) 
(93) 



For short flanges the first term in (90) must be re- 
tained. 

The results expressed by equations (85) to (93) for 
flanges with a simply supi)orted edge at ?/ = 0 could 
have been obtained directly through the approximate 
method by assuming, for cxami)le, instead of (74) that 

W = Ay. 

This was done for isotropic flanges by Timoshenko. 
(Reference 21, p. 396.) 

Effect of elastic giving of material at the base of the 
flange. 

In obtaining the preceding results the lack of fixity 
of the edge y = 0 was ascribed to the twisting of the 



J. 
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column, either as a whole or in segments, in conse- 
quence of the moments applied at this edge by the 
deformation of the flange. Actually, however, the 
material at the base of the flange yields elastically 
under the action of these moments so that the angle 
of rotation of the section is less than {dw/dy)y,= oy 
the inclination of the flange at its base. Accordingly 
equation (72) should be replaced by 

where r; is some proper fraction. The effect is to 
reduce the strain energy (equation (73)) resulting 
from the twisting of the column. To the reduced Vi 
should be added the energy of deformation of the ma- 
terial at the base of the flange. This portion of the 
energy is relatively small. The result is that Fi, equa- 
tion (73), which was added to V, equation (58), to ex- 
press the whole energy of deformation of the flange 
and column in so far as it arises from the load on the 
flanges, should be reduced. This is equivalent to say- 
ing that e as calculated by (79) from the torsional ri- 
gidity of the section should be reduced. 

For flanges of wood in which the grain is longitudinal, 
such reduction in the fixity coefficient is very great. 
This is due to the extremely small relative value of the 
modulus of elasticity E in the direction parallel to the 
faces of the flange and perpendicular to its length, 
which ranges from Yn of the modulus in the longitudinal 
direction in quarter-sawn flanges of spruce to Koo of this 
modulus in flanges in which the gi'owth rings make an 
angle of 45'' with the faces. The tests show^ that, for 
calculated coefficients of fixity of the order of magni- 
tude of 2 and above, the critical stress corresponds to an 
actual fixity of about 0.01. Corresponding reductions 
in the smaller calculated fixity coefficients are ob- 
served but the law that the reduction follows has not 
been determined. 

The practical result of the reduction in fixity because 
of elastic giving is that the condition of a simply 
supported edge at the base of the flange is closely 
approximated when the calculated fixity coefficient is 
small. The material is unable to transmit the bending 
moment from the base of the flange to the body of the 
column, with the result that the flange itself is inclined 
nearly as if it were merely hinged or simply supported 
at its base and consequently a condition in which 
formula (87) is applicable is approached. This 
situation wiU be discussed further in connection with 
the study of two flanged columns with the aid of the 
curves of Figure 23. 

A similar but probably not so great a reduction 
occurs in the calculated fixity coefficients of the 



flanges of structural steel colunms in consequence of 
the elastic giving of the material at the bases of the 
flanges. Practically no data are available for use in 
determining the extent of this reduction. 

Examples of the determination of the critical stresses, 
neglecting the effect of elastic giving at the bases of 
the flanges. 

In the following paragraphs will be explained the 
procedure to follow in applying the results of the 
preceding mathematical analysis, using the fixity 
coefficient as calculated from the torsional rigidity 
of the section and the dunensions of the flange and 
neglecting the reduction in this coefficient that should 
be made to allow for the elastic giving of the material 
at the base of the flange. The method can then be 
applied when the reduced coefficients are known by 
substituting in each case for the fixity coefficient e the 
reduced fixity coefficient e'. 

The method will be first appHed to a column of 
spruce similar to many of those used in the tests. The 
dimensions are shown in Figure 25. The growth 
rings in the single outstanding flange will be assumed 
to make an angle of 45° with the faces of the flange. 
The fixity coefficient is given by 



G 6 



<== 12(1 — (Txy(^yx)T^^^ K 



If 



G_ 1 
E, 15.6 

it follows that 

AVith the given dimensions 
where a = c/b. 



(94) 



6= 1066 -r,- 

a 




Figure 25.— The cross 
section of a wood test 
column with a single 
thin, outstanding flange 
the growth rings of 
which make an angle 
of 45° with the faces 



It is important to observe that 
the coefficient e depends upon 
the half wave length c. This co- 
efficient was computed for a series 
of possible half wave lengths, the 
length of the column being 40 
inches, and the quantity A*, to w^hich 
the corresponding critical stress is 
proportional, was then taken by extrapolation from the 
curves of Figure 23. The results are shown in Table 
XV. The numbers in the last column of the table are 
really estimated, since the values of e concerned are far 
beyond the limits plotted on the curves of Figure 23. 
Through inspection of this column and the curves in 
Figure 23, however, it becomes clear that the flange will 
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break into five half wave lengths, the critical stress be- 
ing 0.12 EJi^lb'^j corresponding to the value 5.33 of the 
ratio cjh. These values agree well with those calcu- 
lated for e= oo. Indeed it is apparent from the be- 
havior of the portions of the curves shown that the 
ordinates rapidly approach their limiting values as e 
increases. 

The approximate method used in calculating the 
curves of Figure 23 gives values of k that are slightly 
too large for the higher values of e. High values of e, 
however, do not occur in cases of practical interest, as 
will shortly be seen. The approximate method may 
therefore be considered entirely satisfactory. 

TABLE XV 

VARIATION OF CRITICAL STRESS WITH NUMBER 
OF HALF WAVE LENGTHS FOR THE 40-INCH COL- 
UMN OF FIGURE 25 



Number of 
half wave 
lengths 


c 


a = c/b 


e 


k 

(estimated) 


8 


5.00 


3. 33 


96.1 


0.18 


7 


5.71 


3. 81 


73.4 


.13 


6 


6. 67 


4. 45 


53.8 


.12 


5 


8.00 


5. 33 


37.5 


. 12 


4 


10.00 


6. 67 


24.0 


.12 


3 


13. 33 


8.89 


13.5 


. 14 


2 


20.00 


13. 33 


6.0 


.20 


1 


40.00 


26. 67 


1.5 


.40 



Table XV was calculated on the assumption that the 
effects of the elastic giving of the material at the base 
of the flange could be neglected. This table indi- 
cated a minimum critical stress of 0.12 Exh^/b^j corre- 
sponding to the value 5.33 of the ratio c/b. Actual 
tests, however, show that the flange wrinkles at a 
stress of 0.07 E:,hyb\ (Part I, equation (5), p. 9.) 
This reduction in the critical stress should be attrib- 
uted to the elastic giving of the material at the base 
of the flange. The curves show that this minimum 
critical stress should be attributed to a fixity coeffi- 
cient in the vicinity of 0.01 and a ratio of c/b of about 7. 
This example is very informing, since it indicates a 
reduction in the fixity coefficient from a number of 
the order of 20.0 to one of the order of 0.01. 

In the example just considered there was onl}^ one 
outstanding flange. If there are N flanges, the fixity 
coefficient as calculated should be divided by A^. 

Consider now the section of column T-25, Figure 5. 
The length of the column is taken as 120 inches. 
The growth rings of the wood will be assumed to 
make an angle of 45° with the faces of the flanges. 
In accordance with equation (79) 



7.58 K b' 

" 4c bh' 



^~ ' Lfi3 ^2~2.736 2 



(95) 



where a = c/b. Proceeding as before Table XVI was 
constructed with the aid of the curves of Figure 23. 



TABLE XVI 

VARIATION OF CRITICAL STRESS WITH NUMBER OF 
HALF WAVE LENGTHS FOR COLUMN T-25 OF 
FIGURE 5 



Number of 










half wave 


c 


a=c/6 




k 


lengths 








8 


If). 00 


4. 67 


0. 1255 


0. 115 


7 


17. 14 


5. 34 


.0960 


.112 


6 


20.00 


6.23 


. 0705 


. 107 


5 


24.00 


7.47 


. 0490 


.105 


4 


30.00 


9. 34 


. 0314 


.103 


3 


40.00 


12, 46 


.0176 


.102 


2 


fiO.OO 


18.68 


.0078 


.107 


1 


120.00 


37. 37 


.0020 


M20 



1 Estimated. 

The values of k in Table XVI indicate that at a 
critical stress of 0.102 EJi^b^ each flange will break 
into three half wave lengths corresponding to a fixity 
coefficient of 0.0176. The tests showed that each 
flange broke into a single half wave length and the 
column twisted at a critical stress of about 0.044 
EJi^l¥, the critical stress for a simply supported edge. 
This means that the calculated fixity coefficient has 
been reduced nearly to zero by the elastic giving of the 
material at the bases of the flanges. 

Failure through twisting or wrinkling. 

When, as in the example just given, the least critical 
stress is associated with a half wave length equal to the 
length of the column, the column fails by twisting 
about its axis. At the base of each flange, as a result 
of the beginning of failure, a torque that is in the same 
sense for the entire length of the column is applied to 
the column as a whole. If a flange breaks into several 
half wave lengths, however, the torques at its base are 
in opposite senses in adjoining half wave lengths and 
consequently oppose one another. 

Practical rules for determining the critical stress, al- 
lowance being made for elastic giving of the material 
at the bases of the flanges. 

In a cruciform section having equal arms and no 
fiUets it appears from equation (95) that a change in 
the dimensions, 6, the outstanding width, and A, the 
thickness of the flange, will not greatly alter the cal- 
culated fixity coefficient e, since if, the torsion constant 
is nearly proportional to b and to K\ (Part I, p. 7.) 

Much the same situation exists in other sections, 
such as L, U, Z, and T, made up of component rect- 
angles, all parts being of equal thickness and having 
no fillets. It appears from the data at our disposal 
that the flanges of such sections may be treated as 
having their bases simply supported. The critical 
stress for long columns of spruce of such sections may 
then be taken as 0.044 E^ h'^lb'\ provided that this 
stress is less than the one that would cause primary 
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failure. If fillets are added to any of these sections 
or if the thickness of the back of a channel is increased, 
for example, the critical stress will increase. Tho 
exact amount of this increase can not be stated, since 
the law by which the calculated fixity coeflicient is 
reduced through the giving of the material at the 
bases of the flanges is not known. Tests indicate, 
however, that the critical stress is increased approxi- 
mately in the ratio of the torsional rigidity of the 
changed section to that of the original section. This 
relation may be taken to hold for spruce until the Hmit- 
ing critical stress 0.07 EJi^b^ is attained. From this 
point as the torsional rigidity increases the critical stress 
remains unchanged. 

As the critical stress increases with increasing co- 
efficient of fixity at the base of the flange, the type of 
failure changes from one through twisting to one 
through wrinkling. The distinction between these 
externally different types of failure does not appear to 
be important, since the one goes over gradually to the 
other. 

For flanges on short columns the critical stresses 
will be higher than those for the long columns just 
considered. 

As previously stated, the foregoing discussion applies 
to flanges of spruce in which the growth rings make an 
angle of 45° with the faces of the flange. Flanges of 
steel or other isotropic material can be treated in a 
similar way through the use of Table XIV and the 
curves of Figure 24. Sufficient experimental data for 
steel columns, to enable the authors to estimate the 
effect of the reduction in the calculated coefficient of 
fixity, have not been found in the literature. 

CONCLUSIONS 

1. Under a compressive load, the critical stress for a 
moderately long flange of spruce, perfectly fixed along 
its base and of thickness h and width 6, is given by 

p = 0.228 E'zJ^ 

when the growth rings are perpendicular to the faces 
of the flange (fig. 21), by 

p = 0.117 E,y, 

when the rings make an angle of 45° with the faces 
(fig. 22), and by 

p = 0.164 E^^2 

when the rings are parallel to the faces. In these 
I'ormulas is Young's modulus in the direction of the 



grain of the wood, which is taken as the direction of the 
length of the flange. 

For a flange of steel the base of which is perfectly 
fixed the critical stress is given by 

V = 1.1^E^, 

when Poisson's ratio is taken as 0..3. 

2. If the base of the flange is simply supported the 
corresponding critical stresses are 

= 0.053 E'xJJ 

and 

2? = 0.044 £',p 

and 

2? = 0.037 E'xp 
for a flange of spruce and 

2^ = 0.385 E'p 

for a flange of steel. Such a condition at the base of 
the flanges is found, for example, in the case of columns 
of L, U, Z, T, and + sections without fillets and 
having parts of the same thickness. 

3. The condition of perfect fixity is not realized in 
practice because of the elastic giving of the material 
at the base of the flange. Tests indicate that the 
upper lunit of the critical stress for moderately long 
flanges of spruce is given by 

2? = 0.07£;.p- 

This is an average value from tests of specimens in 
which the growth rings were at various inclinations 
to the faces of the flanges. For strictly quarter-sawn 
flanges the critical stress would be somewhat higher 
and for plain-sawn ones somewhat lower. The re- 
duction from the values given for flanges with per- 
fectly flxed edges should be attributed to the elastic 
giving of the material at the bases of the flanges. 

Because of the same elastic giving the flxity of 
flanges with partially flxed bases is greatly reduced. 
For such flanges the critical stress ranges from 

p = 0.044£;,J' 

to the upper limit 

2? = 0.07£;x^2- 

Both limiting stresses can be increased somewhat for 
strictly quarter-sawn flanges and should be reduced 
somewhat for plain-sawn ones. 
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Tests on steel flanges were not made. As a result 
of the elastic giving of the material at the base of the 
flange, however, it is probable that the upper limit of 
the critical stress will be found to be considerably less 
than that calculated for a flange with a perfectly 
fixed edge. 

4. The critical stresses for short flanges are greater 
than those given by the preceding formulas. 

5. The critical stresses obtained through use of 
these formulas will be of interest only if they are less 



than those that would cause a primary failure of the 
column under consideration. 



Forest Products Laboratory, 
Forest Service, United States 
Department of Agriculture, 
Madison, Wis., October 15, 1930. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 
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yawing 


pitch 

yaw 









Absolute coefficients of moment 



qhS 



Cm 



gcS 



Cn 



N 
qbS 

4. PROPELLER SYMBOLS 



Angle of set of control surface (relative to neu- 
tral position), 8. (Indicate surface by proper 
subscript.) 



D, Diameter. 

Pj Geometric pitch. 

pIDy Pitch ratio. 

V\ Inflow velocity. 

Vsf Slipstream velocity. 

Ty Thrust, absolute coefficient 
Qy Torque, absolute coefficient Cq- 

5. 

1 hp = 76.04 kg/m/s=-550nb./ft./sec. 
1 kg/m/s = 0.01315 hp 
1 mi./hr.= 0.44704 m/s 
1 m/s = 2.23693 mi./hr. 



P, Power, absolute coefficient Cp==^^3^* 

Cs, Speed power coefficient = -^^2- 
77, Efficiency. 

71, Revolutions per second, r. p. s. 

Effective helix angle = tan"^ (27^} 



Q 

NUMERICAL RELATIONS 



1 lb. = 0.4535924277 kg 
1 kg = 2.2046224 lb. 
1 mi. = 1609.35 m = 5280 ft. 
1 m = 3.2808333 ft. 



